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ABSTRACT 


The  phenomenon  of  cavity  flow  oscillation  is  investigated  to  determine  the 
conditions  for  onset  of  periodic  oscillations  and  to  understemd  the  relationship 
between  the  state  of  the  shear  layer  and  the  cavity  drag.  Experiments  have 
been  performed  in  a  water  tunnel  using  a  4"  axisymmetric  cavity  model  instru¬ 
mented  with  a  strip  heater  on  the  nose  cone  and  pressure  taps  in  and  around 
the  cavity.  A  complete  set  of  measurements  of  oscillation  phase,  amplitude 
amplification  along  the  flow  direction,  distribution  of  shear  stress  eind  other 
momentum  flux  is  obtained  by  meams  of  a  laser  Doppler  velocimeter.  Drag 
measurements  were  made  by  integrating  the  mean  pressure  over  the  solid  sur¬ 
faces  of  the  cavity.  Results  indicated  exponential  cavity  drag  dependence  on  the 
length  of  the  cavity.  A  jump  in  the  cavity  drag  coefficient  is  observed  as  the  cav¬ 
ity  flow  shows  a  bluff  body  wake  type  behavior.  An  independent  estimate  of  the 
drag,  which  is  obtained  by  integration  of  shear  and  mean  momentum  transfer 
terms  over  the  peripheral  area  of  the  cavity,  confirms  the  exponential  depen¬ 
dence  of  drag  on  the  length  of  the  cavity.  Results,  also  reveal  that  the  drag  of 
the  cavity  in  the  non-oscillating  mode  is  less  than  the  case  if  the  cavity  were 
replaced  by  a  solid  surface.  Natural  and  forced  oscillations  of  the  cavity  shear 
layer  spanning  the  gap  are  studied.  The  forced  oscillations  are  introduced  by  a 
sinusoidally  heated  thin-film  strip  which  excites  the  Tollmein-Schlichting  waves 
in  the  boundary  layer  upstream  of  the  gap.  For  a  sufficiently  large  gap,  self- 
sustained  periodic  oscillations  are  observed,while  for  smaller  gaps,  which  do  not 
oscillate  naturally,  periodic  oscillations  can  be  obtedned  by  external  forcing 
through  the  strip  heater.  In  the  latter  case  resonance  is  observed  whenever  the 

forcing  frequency  satisfies  phase  criterion  and  amplitude  exceeds 


certain  threshold  levels,  but  the  phenomenon  is  non-self-supporting.  The  drag 
of  the  cavity  can  be  increased  by  one  order  of  magnitude  in  the  non-oscillating 
case  through  external  forcing.  For  naturally  occurring  oscillations,  it  is  possible 
for  two  waves  to  co-exist  in  the  shear  layer  (natural  and  forced).  Also,  it  is  pos¬ 
sible  to  completely  eliminate  mode  switching  by  applying  external  forcing.  For 
the  first  time  a  test  is  performed  to  ceuicel  or  dampen  the  amplitude  of  Kelvin- 
Helmholtz  wave  in  the  cavity  shear  layer.  This  is  done  through  introducing  an 
external  perturbation  with  the  same  frequency  of  the  natural  component  but 
having  a  different  phase.  Reduction  by  a  factor  of  2  is  obtained  in  the  amplitude 
of  the  oscillation. 
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quency  at  four  different  forcing  power  levels;  Case  A 
6/9o  =  66,  U,  =  23cm Aec. 

Variation  of  the  maximum  shear  stress  in  the  shear 
layer  with  forcing  frequency,  constant  power. 
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Distribution  of  the  maximum  shear  stress  srith  the 
streamwise  distance  along  the  cavity  shear  layer;  a) 
three  different  forcing  frequencies,  constEmt  power; 
b)  2  different  forcing  powers,  constant  frequency; 
6/(9b  =  77. 

Flow  visualization  of  a  naturailly  non*oscillating  cavity 
flow,  b/Ba  =  68: 

(a)  unforced 

(b)  forced  at  F  =  4  Hz,  simulating  Mode  I  of  oscillations. 


Oscilloscope  traces  of  forced  fluctuations. 

(a)  turn  on:  (b)  turn  off; 

b/Bfi  =  86.F  =  7Hz,  2. 

Distribution  of  (a)  term,  (b)  shear  stress  along 

the  y  s  0  line  in  the  non-oscillating  mode; 

O unforced,  4  forced  at  6  Hz. 

Velocity  fluctuation  spectra,  presenting  interaction  of 
the  natural  and  forced  frequencies  in  the  cavity 
shear  layer  for  different  power  levels  (P^),  (h)  scope 
trace  of  Case  (d). 

Possible  interactions  with  the  cavity  flow  through 
external  forcing. 

1}  Mapping  of  Mode  II  into  Regions  1  and  2. 

2)  Mapping  of  Mode  III  into  Region  3. 

3)  Extension  of 

4)  Delay  of  the  wake  Mode. 

5)  Simulation  of  Mode  I. 

Schematic  of  the  Cancellation  Set-up 

Damping  of  natural  oscillations  through 
external  forcing. 

(a)  natural  oscillation 

(b)  strip  heater  on. 
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Symbol  Description 

A  frontal  area  for  the  model 

A I  cavities  peripheral  area 

Ag  downstream  face  area  (=A) 

As  stings  peripheral  area 

A4  upstream  face  area 

b  cavity  width 

minimum  cavity  width  for  self-sustained  oscillations 
b„  cavity  width  for  transition  to  the  wake  mode 
c  speed  of  sound 

Co  cavity  drag  coefficient 

Cos  sting  drag  coefficient 

Cop  pressure  drag  coefficient  of  the  cavity 

drag  coefficient  based  on  uv  term 
Cinr  drag  coefficient  based  on  edge  velocity  d3m£unic  pressure 
Cp  pressure  coefficient  based  on  edge  velocity  dynamic  pressure 

^^max  maximum  locaJ  pressure  coefficient  of  the  cavity 

d  model  diameter 

dr  length  element  zdong  the  model  radius 

dx  length  element  along  the  cavity  width 

ds  sting  diameter 

F  frequency  of  the  cavity  flow  oscillation 

h  two-dimensional  cavity  depth 

H  boundary  layer  shape  factor 


K  LDV  conversion  factor 

Kg  LDV  conversion  factor,  blue  channel 

Kff  LDV  conversion  factor,  green  channel 

2L  direction  cosine  for  the  cavity  surface 

P  static  pressure 

Pg  free  stream  static  pressure  oustide  the  upstream  boundary  layer 

r  radius  from  center  line  of  the  model 

R  model  radius 

R,  Reynolds  number  based  on  model  diameter 

Re,.  Reynolds  number  based  on  boundary  layer  displacement  thickness 

Re®  Reynolds  number  based  on  boundary  layer  momentum  thickness 

S  stream”~ise  distance  over  the  model  nose 

SH\  location  of  first  strip  heater 

SH2  location  of  second  strip  heater 

SHz  location  of  third  strip  heater 

u  stream’.vise  velocity  component 

C4  convection  or  phase  velocity  of  the  perturbation 

Ug  free  stream  velocity  outside  the  upstream  boundeu’y  layer 

V  transverse  velocity  component 

w  strip  heater  width 

X  streamwise  direction 

y  outward  direction  of  the  model  perpendicular  to  x-z  plane 

i/,a  locus  of  -^=.5 

t/g 

z  vertical  direction  of  the  tunnel 

a  angle  between  laser  beams 

d'  boundary  layer  displacement  thickness 

phase  difference  between  two  signals 


time  between  data  points  in  LDV  signal 

time  shift  in  cross  correlation  calculation 

boundary  or  shear  iayer  momentum  thickness 

wave  length  of  perturbations 

acoustic  wave  length 

fluid  density 

viscous  shear  stress 

overall  phase  difference  between  two  comers 

Doppler  frequency 

upstream  condition 

minimum  value 

maximum  value 

time  average  or  mean  value 

fluctuation  about  the  mean 

rms  value 

vector  quantity 

free  stream  condition  outside  the  upstream  boundary  layer 


Chapter  1 


Introduction 


1.1.  Background 

Cavity  flow  phenomena  occur  in  a  veu'iety  of  structures;  slots  between  mov¬ 
able  parts  of  control  surfaces  of  aircraft,  ships,  tandem  eirrangements  of  bluff 
bodies  such  as  adjacent  tall  buildings  or  tractor-trailer  combinations,  gas 
d3niamic  lasers,  hydraulic  gates,  and  control  valves.  The  presence  of  cavities  can 
give  rise  to  drastic  changes  in  the  steady  and  unsteady  loading  character  of 
these  structures.  The  flow  past  cavity-type  geometries  (cut-outs)  may  introduce 
large  amplitude  velocity  and  pressure  fluctuations  in  and  around  the  cavity 
which  are  often  damaiging  to  the  structures.  Flow  oscillation  over  gaps  between 
insulation  tiles  used  on  the  space  shuttle  outer  surface,  caused  by  evaporation 
of  the  flller  material  at  reentry,  is  suspected  of  being  responsible  for  the  loss  of 
a  portion  of  the  tiles.  Emission  of  acoustical  noise  is  another  unwanted  aspect 
of  cavity  flow-type  oscillation.  On  the  other  hand,  flow  over  cavities  can  be 
effectively  utilized  to  alter  heat  transfer  (White  1971)  and  drag  characteristics  of 
aerodynamic  configurations  (Koenig  1976). 

Charwat,  et  al.  (1961)  categorized  cavity  flows  into  open  or  closed  types: 
"open"  when  the  separated  shear  layer  spans  the  gap, while  "closed"  refers  to 
flows  where  the  separated  shear  layer  reattaches  on  the  bottom  of  the  cavity 
and  separates  again  ahead  of  the  downstream  corner.  For  supersonic  speeds 
and  for  subsonic  speeds  with  turbulent  boundary  layers,  closed  cavities  are 
found  to  have  a  width-to-depth  b/h  ratio  of  more  than  11  (Charwart,  et  al. 
1961).  Sarohia  (1975)  shows  a  demarcation  line  between  an  open  and  closed 
cavity  to  be  around  b/h  v  7-8  for  low  subsonic  speeds  with  laminar  upstream 


boundary  layers.  Furthermore,  open  cavity  flows  are  divided  into  shallow  and 
deep  cavities,  depending  on  the  t5rpe  of  acoustic  resonance  present  in  the  flow. 
When  an  acoustic  resonance  is  present  in  the  flow,  the  standing  acoustic  wave  is 
oriented  longitudinally  for  shallow  cavities  (6/fc  >  1).  The  standing  wave  is  of 
transverse  type  for  deep  cavities  (b/h  <  1)  (Charwat  et  al..  1961.  Plumblee  et  al.. 
1962.  East  1966.  Tam  and  Block  1978).  Such  a  definition  for  the  shallow  and 
deep  cavities  loses  its  meaning  in  the  absence  of  the  acoustic  resonance  in  the 
cavity.  Tani  et  al..  (1961)  define  cavity  flows  as  deep  and  shallow  cavities  based 
on  their  respective  low  amd  high  pressure  drag  values.  The  demarcation  line 
based  on  Tani's  definition  is  about  fa  A  »  1.4. 

Sketch  1  illustrates  the  physically  significant  parameters  of  the  cavity  flow. 


where  U,  is  the  upstream  edge  velocity  outside  the  boundary  layer,  6*0  is  the 
displacement  thickness  and  Sq  is  momentum  thickness.  For  fixed  flow  parame¬ 
ters  ( £/,.  <5*o)'  cavity  flow  shows  a  strong  oscillation  of  the  shear  layer  when 
the  width  b  exceeds  a  critical  minimum  width,  6min-  (Karamcheti  1956, 
Sarohia  1975,  Rockwell  aind  Naudascher  1978).  The  existence  of  this  minimum 


width  for  the  onset  of  oscillation  is  linked  to  the  nature  of  the  instability  growth 
in  the  shear  layer.  Previous  investigators  noticed  that  the  main  driving  force  of 
the  acoustic  resonance  condition  in  the  cavity  was  the  amplification  of  instabili¬ 
ties  by  the  respective  cavity  shear  layer.  In  the  shear  layers,  the  final  grovrth 
and  saturation  of  the  instability  waves  result  in  roll-up  of  the  shear  layer  and 
appearance  of  large  vortical  structures  (Brown  emd  Roshko  1974).  The  velocity 
fluctuation  associated  with  the  cavity  flow  oscillations  can  be  attributed  to  the 
presence  of  these  large  vortical  structures  (Sarohia  1975,  Rockwell  and  Knisely 
1979).  In  his  study  of  the  frequency  selection  process,  Sarohia  (1975)  concluded 
that  the  frequency  at  which  the  cavity  chooses  to  oscillate  is  such  that: 

(1)  certain  phase  criteria  necessary  for  the  self-sustainment  of  the  flow  are 

satisfied;  and 

(2)  maximum  integrated  amplification  through  the  cavity  shear  layer  occurs. 

The  cavity  flow  oscillation  was  described  as  a  self-sustained  type.  Rockwell  and 
Naudascher  (1970)  suggest  the  following  chain  of  events  to  be  responsible  for 
the  self-sustaining  nature  of  the  oscillation:  impingement  of  organized  vortical 
structures  upon  the  downstream  corner:  resultant  upstream  influence  (inter¬ 
preted  as  Biot-Savart  induction);  reception  and  conversion  of  this  upstream 
influence  into  velocity  perturbation  by  the  sensitive  region  of  the  shear  layer 
close  to  the  upstream  corner:  amplification  of  these  fluctuations  by  the  shear 
layer;  growth  into  vortex  structures.  Consideration  of  the  acoustic  speed 
becomes  important  in  the  nature  of  the  upstream  propagating  waves  for  high 
Mach  number  flows  or  when  the  impingement  length,  b,  is  of  the  same  order  of 
magnitude  as  the  acoustic  wave  length  (b  »  \a)  (Heller  and  Bliss  1975,  Rossiter 
1964,  and  Rockwell  and  Naudascher  1978). 


The  self-sustained  nature  of  the  cavity  flow  oscillation  is  similar  to  that  of  the 
jet-edge  oscillation  (edge  tone  phenomenon).  Self-sustained  oscillation  of  cavity 
flow  has  the  following  properties: 

(1)  for  a  fixed  cavity  depth  and  fixed  flow  parameters  ( Ue_0o6^  there  exists  a 
minimum  width  (bmin)  below  which  no  oscillation  occurs; 

(2)  for  fixed  cavity  depth  (h)  and  fixed  boundary  layer  thickness  there  exists 
a  minimum  Mach  number  or  free  stream  velocity  ( Ue)  below  which  no  oscillation 
occurs ; 

(3)  for  b  >  or  Ue  >  Ug^  and  fixed  boundary  layer  thickness  the  fre- 

Ug 

quency  of  oscillation  decreases  as  the  parameter  — —  decreases. 

b 

At  certain  values  of  the  frequency  undergoes  a  jump  to  a  higher  mode 

of  oscillation  with  a  corresponding  higher  frequency.  Hysteresis  is  observed  in 

u 

the  transition  region.  In  the  new  mode,  frequency  also  decreases  as  — ^ 

b 

decreases.  Sarohia  (1975)  has  found  that  for  fixed  values  of  ^4.  b,  and  h,  there 
exists  a  maximum  boundary  layer  thickness  or  do'maxj  above  which  no 

oscillation  occurs. 

The  ladder-type  behavior  of  frequency  versus  — -  (which  has  the  dimension 

h 

of  frequency)  is  remarkably  similar  for  self-sustained  flows  and  externally 
forced  free  shear  layers  (Ho  and  Huang  1982).  However, there  is  a  significant 
difference  between  the  nature  of  perturbations  introduced  at  the  sensitive 
region  of  the  shear  layer  by  the  downstream  corner  and  that  of  an  external 
forcing  source.  External  forcing  subjects  the  shear  layer  to  a  constant  ampli¬ 
tude  perturbation  at  one  frequency  But  flows  excited  by  the  impingement  of 
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the  shear  layer  at  the  corner  are  subjected  to  variation  in  strength  of  events  at 
the  corner  (Knisely  1980).  Often  this  variation  gives  rise  to  the  amplitude  modu¬ 
lation  of  the  fluctuations. which  shows  as  a  low  frequency  component  in  the 
spectrum  of  fluctuations.  Knisely  (1980)  shows  that  transverse  positioning  of 
the  successive  vortices  as  they  approach  the  downstream  corner  yields  vary¬ 
ing  incident  structure  severing  degrees.  It  is  not  clear  what  sort  of  mechan¬ 
ism  is  responsible  for  the  presence  of  these  low  frequencies  in  the  cavity  shear 
layer.  Also,  the  precise  role  of  the  downstream  corner  as  a  source  of  external 
forcing  for  the  shear  layer  and  generation  of  different  frequency  components  of 
oscillation  is  not  fully  understood. 


As  mentioned  before,  the  presence  of  cavities  (gaps)  can  be  used  effectively  to 
change  the  drag  characteristics  of  bluff  bodies.  A  novel  experiment  by  Roshko 
and  Koenig  (1978)  demonstrated  that  by  properly  choosing  dimensions,  a  gap 
generated  by  placing  a  disk  in  front  of  a  circular  cylinder  can  drastically  reduce 
the  drag  coefficient  of  the  system.  A  similar  observation  for  different 
configurations  of  bluff  bodies  in  tandem  was  made  by  Mair  (1965)  and  Zdravko- 
vich  and  Pridden(1977) .  It  is  unclear  into  which  category  these  oscilla¬ 
tions  fall.  Scattered  information  about  the  frequency  of  the  shear  layer  oscilla¬ 
tion  shows  behavior  characteristic  of  self-sustained  flow  for  the  low  drag  cases 
of  flow  over  gaps  generated  by  bluff  bodies  in  tandem. 


1.2.  Statement  of  the  Problem 

A  large  portion  of  cavity  flow  investigations  are  devoted  to  the  study  of  insta¬ 
bility  aspects  of  this  problem.  Quite  independently,  a  second  group  of  studies 
can  be  found  which  deal  only  ynth  the  pressure  field  and  drag  characteristics  of 
this  flow  (Roshko  1955,  Tani  1961,  Fox  1964,  Koenig  1978).  To  our  knowledge  no 
investigator  has  tried  to  fill  the  gap  between  the  first  and  the  second  group  of 
experiments.  It  was  this  gap  that  motivated  the  present  investigation  into  the 
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role  of  instabilities  in  drag  and  pressvire  field  formation  of  cavities. 

Another  important  aspect  of  cavity  fiow  study  is  to  develop  a  technique  to 
control  the  frequency  and  amplitude  of  oscillations.  Previous  work  involving 
control  of  the  cavity  flow  oscillation  has  dealt  Yhth  changing  the  geometry  of  the 
downstream  or  upstream  corner  of  the  cavity  by  rounding  the  corners  (Ethem- 
babaoglu  1973)  or  by  placement  of  external  objects  in  the  cavity  such  as  flaps, 
deflectors  or  perforated  plates  (Heller  and  Bliss  1975  and  Ball  1957).  Usually 
these  devices  introduce  their  own  complexities  into  the  problem.  There  was  a 
need  for  a  clean  technique  that  can  control  the  oscillations  by  purely  introduc¬ 
ing  external  perturbation  to  the  sheau"  layer. 

The  following  were  objectives  set  for  the  present  investigation. 

(1)  To  understand  the  relation  between  the  state  of  the  shear  layer  and  the 
drag  of  the  cavity. 

(2)  To  determine  the  influence  of  initial  conditions  on  the  onset  of  self- 
sustained  oscillation  and  to  obtain  a  better  understanding  of  the  feed-back 
mechanism. 

(3)  To  develop  and  use  a  technique  which  is  capable  of  controlling  (enhancing 
or  damping)  established  oscillations. 

1.3.  Approach 

The  connection  between  the  state  of  the  sheeu"  layer  and  the  drag  of  the  cav¬ 
ity  system  can  be  made  by  studying  the  quantities  that  are  responsible  for  the 
transfer  of  the  momentum  in  and  out  of  the  cavity  {uW  andiliJ).  These  quanti¬ 
ties  carry  the  necessary  information  about  the  nature  of  the  instabilities  in  the 
shear  layer  and  play  an  inherent  role  in  the  formation  of  the  pressure  distribu¬ 
tion  inside  the  cavity  system.  Thus,  accurate  measurements  of  these  quantities 
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for  different  width-to-height  ratios  (b/h)  in  conjunction  irith  the  pressure  dis¬ 
tribution  measurements  would  provide  the  necessary  information  for  fulfilling 
the  first  objective.  External  forcing  by  strip  heating  was  used  to  study  the  feed¬ 
back  mechanism  and  control  of  cavity  fiow  oscillation, thus,  fulfilling  the  second 
and  third  objectives. 


Chapter  2 


Experimental  Details 


2.1.  Water  Channd 

The  experiments  were  performed  in  a  newly  constructed  free  surface  water 
tunnel,  which  is  located  at  the  Jet  Propulsion  Laboratory,  California  Institute  of 
Technology.  This  tunnel  is  specially  designed  to  attain  a  low  turbulence  level  for 
the  purpose  of  carrying  out  sensitive  instability-t3T>e  experiments.  The  essential 
features  of  the  free  surface  tunnel  are  shown  in  figure  1.  Water  is  recirculated 
in  a  closed  loop  by  a  50  horsepower  axial  fiow  pump,  located  5m  below  the  free 
surface,  in  the  lower  leg  of  the  tunnel.  The  desired  fiow  speeds  are  achieved  by 
throttling  fiow  through  an  on-line  and  a  by-pass  valve.  The  upper  portion  of  the 
tunnel  consists  of  a  diffuser,  a  straight  wall  section  comprised  of  fiow  manipula¬ 
tors  and  a  test  section.  The  high  speed  flow  from  the  pump  is  then  diffused  by 
the  22°  half-angle  diffuser.  In  order  to  prevent  severe  flow  sepeiration,  the  effect 
of  an  unfavorable  pressure  gradient  is  reduced  by  matched  pressure  drops  of 
three  perforated  plates  installed  in  the  diffuser  section.  Before  entering  the  test 
section,  the  flow  passes  through  different  flow  manipulators  and  a  final  4:1  fifth- 
order  polynomial  contraction.  Each  manipulator  and  contraction  is  designed  to 
break  up  and  suppress  different  scales  of  the  turbulence. 

The  test  section  is  45.7  cm  (16  inches)  wide  and  63.5  cm  (25  inches)  deep. 
The  wetted  depth  can  be  increased  up  to  56  cm  (22  inches).  The  test  section  is 
bounded  by  glass  windows  on  sides,  bottom  and  top  to  permit  visual  and  optical 
study  of  the  flow  throughout  its  1.63  m  (6  feet)  length. 


The  existing  arrangement  made  it  possible  to  achieve  a  maximum  velocity  of 
1  m/sec  in  the  test  section.  Turbulence  levels  and  velocity  profiles  at  20  cm 
from  the  entrance  to  the  test  section  and  outside  of  the  wall’s  boundary  layer 
were  obtained  using  hot  film  probes  and  am  LDV  system.  The  uniformity  of  the 
flow  within  a  square  area  (30  cm  x  30  cm)  centered  at  the  center  of  the  test  sec¬ 
tion  was  below  1%  of  the  center  line  mean  velocity.  The  free  stream  turbulence 
level  at  25  cm/sec  on  the  centerline  was  near  .05%.  The  power  spectrum  of 
free  stream  velocity  fluctuations,  a  sample  of  which  is  shown  in  figure  2a,  indi¬ 
cates  most  of  the  energy  was  concentrated  at  frequencies  below  .5  Hz,  which  are 
primarily  below  the  range  of  frequencies  of  interest  in  this  study,  2  to  10  Hz 
(figure  2). 

2.2.  Model 

An  axisymmetric  model,  constructed  of  anodized  aluminum  and  Plexiglas, 
was  used  throughout  the  investigation.  The  model  was  composed  of  a  rear  por¬ 
tion  and  a  front  portion  which  were  connected  by  a  sting  as  shown  schematically 
in  figure  3.  The  rear  portion,  a  hollow  cylinder  10.16  cm  (4  inches)  in  diameter 
and  91.4  cm  (36  inches)  in  length,  was  connected  by  a  strut  to  a  traversing 
mechanism  which  was  installed  on  top  of  the  test  section.  This  allowed  for 
movement  of  the  model  up  and  down  and  also  along  the  test  section  axis.  The 
front  portion  consisted  of  a  half  portion  2:1  ellipsoidal  cone  with  a  10.16  cm  (4 
inch)  minor  diameter.  The  ellipsoidal  cone  was  extended  in  the  downstream 
direction  by  a  straight  section  which  edso  was  10.16  cm  (4  inches)  in  diameter 
and  5.06  cm  (2  inches)  in  length. 

The  sting  was  made  from  a  1.27  cm  (1/2  inch)  O.D,  thick  wall  stainless  steel 
tube.  Twelve  brass  tubes  were  installed  inside  of  its  walls  to  carry  pressure  sig¬ 
nals  from  the  front  nose  while  the  interior  of  the  sting  was  devoted  to  dye  tubes 
and  electrical  wires  for  the  strip  heaters.  The  front  end  of  the  sting  was  fixed  to 


the  inside  wedl  of  the  front  nose  while  the  rear  end  was  connected  to  a  vertical 
rod  by  means  of  a  90“  geeiring  mechzmlsm.  The  vertical  rod  was  then  brought 
out  of  the  model  through  the  strut.  By  manual  rotation  of  the  vertical  rod,  the 
sting  could  be  moved  longitudinally  relative  to  the  rear  portion  of  the  model. 
With  the  aid  of  this  mechanism,  plus  a  dial  and  a  helipot  at  the  end  of  the  verti¬ 
cal  rod,  it  was  possible  to  change  and  read  the  cavity  length.  The  solid  inner 
faces  of  the  cavity  were  equipped  with  14  pressure  taps  on  the  downstream  face 
and  8  pressure  taps  on  the  upstream  face.  The  diameter  of  each  tap  opening 
was  .1  cm  (.041  inch  ). 

2.3.  Pressure  Measurements 

The  low  velocities  in  the  range  of  interest  introduced  the  well  known  difficulty 
of  measuring  low  differential  pressures  under  high  head.  A  free  stream  velocity 
of  25  cm/sec  generated  a  dyneunic  head  of  .38  cm  (.15  inch)  of  water,  while  a 
typical  value  of  total  head  at  the  point  where  the  pressure  transducer  was 
located  was  63.5  cm  (25  inches)  of  water.  Further,  since  differential  pressures 
inside  the  cavity  are  much  smaller  than  for  a  free  stream  dynamic  head,  the 
demand  on  the  pressure  measuring  system  for  high  resolution  and  repeatability 
is  very  severe. 

Another  common  problem  in  pressure  measurements  with  water  as  a  medium 
arises  from  the  presence  of  air  bubbles  in  the  system,  which  may  cause  errone¬ 
ous  pressure  readings.  The  following  approach  was  taken  to  solve  the  air  bubble 
problem  and  to  achieve  high  resolution.  Each  line,  Tygon  tube  ID  =  .16  cm  (1/16 
inch),  from  the  model  pressure  taps  or  from  the  pitot-static  probe  was  con¬ 
nected  to  a  common  drain  through  a  toggle  valve  (figure  4).  Each  line  then 
extended  and  manifolded  into  a  24  port  fast  switch  (Scanivalve  t3pe  W2460). 
The  common  end  of  the  Scanivalve  was  connected  to  the  low  pressure  side  of  a 
pressure  transducer  (Validyne  type  Wet-Wet  DP-103-12)  which  was  able  to 


measure  differential  pressures  in  the  range  of  0  to  ,762  cm  (.3  inch)  of  water 
with  a  resolution  of  1%  of  the  total  reuige.  The  high  pressure  side  of  the  trans¬ 
ducer  was  connected  to  the  total  pressure  tube  in  the  test  section,  which  was 
chosen  as  the  reference.  The  common  drain  was  connected  to  a  vacuum  line  and 
was  used  as  a  sink  into  which  air  bubbles  could  rise.  Also,  it  was  used  to  fill  the 
lines  from  the  model  to  the  toggle  valves.  A  special  calibration  tank  was 
designed  to  provide  a  differential  head  of  0  to  10.16  cm  (4  inch)  of  water  with  a 
resolution  of  .00254  cm  (.001  inch)  of  water.  One  port  of  the  Scanivalve  was 
devoted  to  the  low  head  side  of  the  calibration  tank,  while  the  higher  head  side 
was  connected  to  the  high  pressure  side  of  the  transducer  through  a  3-way 
valve.  A  schematic  of  the  pressure  sensing  system  is  given  in  figure  4.  The  aver¬ 
age  height  of  the  calibration  tank  was  set  close  to  that  of  the  free  surface.  This 
was  done  to  prevent  overpressurizing  the  transducer  during  switching.  Table  1 
gives  different  operational  positions  of  the  valves.  The  pressure  measurement 
process  was  automated  by  using  an  HP  data  logger.  A  proper  software  program 
was  written  for  the  data  logger  in  order  to  actuate  the  Scanivalve  and  to  digitize 
the  analog  signal  from  the  transducer.  Prior  to  each  run  the  digitization  rate, 
the  lag  time  between  each  point,  and  the  averaging  time  were  specified  to  the 
machine.  After  each  run,  the  program  stored  the  pressure  distributions,  and 
the  calculated  drag  coefficient  as  well  as  other  parameters  of  each  run  on  a 
floppy  disk. 

2.4.  Frequency,  Phase  and  Correlation  Measurements 

Velocity  fluctuation  measurements  were  made  with  a  TSI  hot -film  probe 
(1277-lOA)  using  a  TSI  1050C  anemometer  in  conjunction  with  a  TSI  true  rms 
meter  (1076).  The  output  of  the  hot  film  anemometer  was  fed  to  a  digital  power 
spectrum  analyzer  (SL20).  A  bandwidth  of  .1  Hz  was  used  for  locating  and  dis¬ 
tinguishing  frequency  content  of  the  velocity  signals  on  line.  A  hard  copy  of  the 
spectrum  was  obtained  by  either  taking  pictures  of  the  spectrum  analyzer 


screen  or  by  plotting  the  spectrum  on  an  x-y  plotter  (HP  model  7001 A). 

For  phase  angle  and  cross  correlation  measurements  a  second  hot-fllm  probe 
was  mounted  on  a  traversing  mechanism  which  was  able  to  traverse  the  entire 
width  of  the  cavity.  Throughout  this  experiment  the  first  probe  was  used  as  a 
reference  probe.  This  probe  was  placed  above  the  upstream  corner  at  x  =  b  and 
I/.95.  the  latter  being  the  y-location  of  UAJ,  =  .95.  Phase  measurements  were 
made  by  moving  the  second  probe  relative  to  the  reference  probe  along  the  edge 
of  the  shear  layer  where  U/U,  =  .95  {y  =yjis)-  The  signals  from  the  reference 
and  the  second  probe  were  fed  to  a  digital  probability  and  correlation  analyzer 
(SAICOR  43A).  Throughout  the  phase  angle  or  cross  correlation  measurements  a 
sampling  rate  of  200  per  second,  limited  by  400  samples  and  a  corresponding 
Nyquist  frequency  of  about  30.  was  used.  Again  a  hard  copy  of  the  output  was 
obtained  on  an  x-y  plotter.  See  Appendix  A  for  the  details  of  phase  angle  meas¬ 
urement. 

2.5.  Velocity  Heosorements  by  Laser  Doppler  Velocimetry 

Due  to  sensitivity  of  the  cavity  flow  oscillations  to  the  presence  of  external 
objects  in  the  flow,  hot-film  anemometry  is  not  suitable  for  velocity  measure¬ 
ments.  Severe  interference  of  the  probes  with  cavity  flow  mode  switching  is 
observed  in  the  present  study.  Existence  of  large  regions  of  reversed  flow  is 
another  reason  for  making  hot-film  anemometry  unappealing  for  this  flow.  The 
only  instrument  currently  available  to  measure  velocity  field  non-intrusively  and 
unambiguously  in  such  a  flow  condition  is  the  laser  Doppler  velocimeter  (LDV). 
Other  useful  properties  of  the  LDV  include  the  small  sampling  volume,  high  fre¬ 
quency  response  and  linearity  of  the  instrument.  In  this  study,  the  velocity  field 
of  the  cavity  flow  was  measured  by  a  two-component,  two-color,  frequency 
shifted  laser  Doppler  velocimeter  (TSI  9100-7  system).  It  is  a  dual  beam  (for¬ 
ward)  scatter  configuration  and  processor,  which  was  operated  in  single-burst 


mode.  The  LDV  processor  operated  in  conjunction  with  a  minicomputer  system 
(PDP-11),  which  was  also  used  for  processing  the  results.  A  forward  scatter 
configuration  was  used  in  these  experiments,  vrith  the  transmitting  and  receiv¬ 
ing  optics  located  at  the  bottom  and  top  side  of  the  facility  as  shown  in  figure  5. 

The  physiced  arrangement  of  the  LDV  optics  relative  to  the  facility  and  infor¬ 
mation  regarding  focal  volume  are  given  in  figure  6  and  table  2.  The  transmit¬ 
ting  and  receiving  optics  are  mounted  on  an  optical  table  with  their  optical  axes 
perpendicular  to  the  flow  direction.  The  focal  volume  has  its  longitudinal  axis  in 
the  z  direction  of  the  facility.  The  streamwise  velocity  component  is  measured 
by  the  fringe  system  formed  in  the  z-x  plane  of  the  tunnel  .while  velocity  com¬ 
ponents  perpendiculeir  to  the  mean  free  stream  are  measured  by  the  fringe  sys¬ 
tem  formed  in  ^he  z-y  plane  of  the  channel.  The  velocity  bias  technique  (Bragg 
cell)  is  used  to  detect  velocity  reversals.  For  details  of  the  optical  setup  and 
Bragg  shifting  see  Appendices  B  and  C,  respectively. 

The  focal  volume  can  be  moved  within  the  test  section  by  means  of  remotely 
operated  stepping  motors.  Motion  on  a  fixed  horizontal  plane  was  obtained  by  a 
corresponding  motion  of  the  optical  table.  Motion  along  the  z-axis  resulted 
from  the  displacement  of  the  focussing  lens  and  receiving  optics  along  their  opt¬ 
ical  axis.  The  table  position  was  controlled  to  a  resolution  better  than  100  /xm. 

2.5.1.  Signal  Processor  and  Data  Acquisition.  The  measurement  of  the  velo¬ 
city  can  be  obtained  only  after  suitable  processing  of  the  outputs  of  the  pho¬ 
tomultiplier  tubes.  Twc  counter  processors  were  used,  one  for  each  velocity 
component.  They  were  operated  in  the  single-burst  mode,  for  which  only  one 
particle  is  present  in  the  focal  volume  when  the  measurement  is  made.  In  addi¬ 
tion,  the  system  was  operated  in  coincident  mode,  in  which  the  output  of  the 
processors  is  tested  for  coincidence  in  time  so  that  the  two  components  of  the 
velocity  vector  as  determined  by  the  optical  set-up  are  measured  on  the  same 


particle. 


The  signal  processing  is  accomplished  in  four  stages:  analog  processing,  burst 
processing,  coincidence  test  and  transfer  to  the  computer.  The  einalog  process¬ 
ing  involves  suitable  filtering  and  amplification  of  the  signal  for  optimum  signal- 
to-noise  ratio.  Only  low  pass  filters  were  used.  A  built-in  filter  in  the  Breigg  cell 
circuitry  removes  pedestal  noise  (low  frequencies).  The  low  pass  filter  is  used 
for  removal  of  high  frequency  noise.  Cut-off  values  from  30  kHz  to  100  kHz  were 
used,  depending  on  fiow  velocity.  The  output  of  the  filter  is  processed  by  an 
amplifier.  In  the  present  experiments  the  gain  was  set  at  2  or  4  depending  on 
the  condition  of  the  signal.  A  detailed  description  of  the  valid  burst  processing 
is  given  in  Appendix  D. 

2.5.2.  Operating  Conditums.  A  t5rpieal  signal-to-noise  ratio,  as  observed  from 
high  sampling  rate  measurements  in  a  quiescent  fluid,  was  in  excess  of  100.  An 
equivalent  background  noise  level  of  .4%  and  .6%  was  detected  for  the  stream- 
wise  and  the  transverse  velocity  fluctuations  .respectively.  When  the  signals 
from  two  channels  were  correlated  for  the  purpose  of  shear  stress  measure- 

ment,  the  corresponding  noise  level  for  "v  /  ^  measurements  was  .07%.  The 

▼  t/, 

better  noise  level  value  was  due  to  the  fact  that  much  of  the  uncorrelated  noise 
from  two  independent  processors  was  eliminated  during  the  correlation  process, 
resulting  in  an  extremely  low  noise  level  for  shear  stress  measurement. 

A  time  window  of  50  fisec  was  used  in  all  velocity  measurements.  The  ampli¬ 
tude  threshold  level  and  the  gaun  of  the  amplifier  were  adjusted  to  maintain  the 
data  rate  between  60  to  200  particles / sec,  which  corresponds  to  a  minimum 
Nyquist  frequency  of  10.  It  was  found  under  the  present  experimental  condi¬ 
tions  that  particles  less  them  15  fim  had  sufficient  response  frequency.  Thus  .for 
the  purpose  of  LDV  measurements,  particles  larger  than  15  ;im  were  filtered 
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from  the  water. 

2.5.3.  Operating  Procedures.  Prior  to  every  experimental  run  with  the  LDV 
s)rstem,  a  variety  of  tests  were  conducted  to  verify  adequate  operation  of  the 
system,  e.g.,  the  transmitting  optics  alignment  and  overlapping  of  the  green  and 
blue  focal  volumes.  The  alignment  of  the  receiving  optics  was  verified  with  the 
facility  in  operation.  Adjustment  of  the  location  of  the  pin-hole  and  photomulti¬ 
plier  assembly  was  often  necessary  for  maximum  data  acquisition  rate  at  the 
operating  conditions  of  the  water  channel. 

After  the  high-pass  and  low-pass  filters  were  set  to  their  predetermined 
values,  the  gain  of  the  amplifier  stcige  was  adjusted  for  optimum  signal-to-noise 
ratio.  After  these  tests  were  completed,  the  actual  measurements  were  ini¬ 
tiated.  A  16,000  word  buffer  was  obtained  at  every  point.  The  data  were  stored 
temporarily  on  disk.  At  the  end  of  the  test  run,  all  the  data  were  transferred  to 
a  magnetic  tape  for  permanent  storage  and  processing  at  a  later  time.  The  size 
of  the  data  buffer  combined  with  the  data  rate  to  give  a  total  sampling  time 
between  15  and  60  sec. This  sampling  time  was  considered  adequate  to  properly 
characterize  mean  turbulent  quantities. 

2.5.4.  Data  Processing.  The  velocity  data  generated  by  the  data  acquisition 
programs  was  processed  using  the  minicomputer  system.  Processing  involved 
several  steps,  e  g.,  error  detection,  computation  of  flow  characteristics  in  the 
frame  of  reference  of  the  optics,  transformation  of  those  quantities  to  the  phy¬ 
sical  frame  of  reference,  and  the  computation  of  integral  quantities  in  the  phy¬ 
sical  reference  frame. 

There  was  one  condition  that  could  have  led  to  a  significant  error  in  the 
measurements.  This  condition  arises  from  transmission  line  noise,  resulting  in 
a  spurious  data  word  being  included  in  the  stream  of  data  as  it  is  transferred  to 


the  computer.  Because  catastrophic  effects  could  result  from  this  condition, 
the  data  processing  programs  had  to  verify  that  the  data  were  in  the  correct 
order.  If  an  error  was  detected,  the  data  point  would  be  discarded.  Less  than 
1%  of  the  data  was  removed  by  this  test.  At  every  measurement  point,  a  fre¬ 
quency  histogram  was  constructed  for  each  channel.  This  histogram  provided 
information  about  the  proper  operation  of  the  system,  i.e.,  filter  setting  and  the 
electronic  noise. 

The  Doppler  frequency  and  the  time  between  data  points  were  utilized  to  cal¬ 
culate  the  fiow  velocity  using  the  appropriate  Doppler  constant.  The  mean, 
second  moment  and  cross-correlation  between  the  two  velocity  components 
were  evaluated.  For  a  detailed  description  of  the  expressions  used  for  the  data 
processing  see  Appendix  E. 

2.6.  Flow  Visualization 

The  model  was  equipped  with  two  dye  injection  tubes  whose  axes  each  made 
an  angle  of  1°  with  local  curvature  of  the  nose  cone  at  the  point  of  the  injection. 
The  injection  points  were  2.54  cm  (1  inch)  off  the  stagnation  point.  Two  visuali¬ 
zation  techniques  were  employed  in  this  study.  In  both  cases,  diluted  coloring 
dye  or  fluorescein  dye  was  gravity  fed  to  the  injection  tubes  and  then  discharged 
into  the  boundary  layer  over  the  cavity. 

In  the  first  method,  using  blue  dye,  a  uniform  background  illumination  was 
provided  specifically  for  photographic  purposes.  In  the  second  technique 
fluorescein  dye  was  used  and  excited  by  a  laser  sheet.  The  laser  sheet  was 
obtained  from  a  2-watt  argon-ion  laser  with  a  thickness  less  than  1mm.  The 
technique  was  used  first  by  Dewey  (1976).  Photographs  of  the  resulting  flow  pat¬ 
terns  were  made  using  a  35  mm  SLR  camera  with  a  focal  length  F3.6  lens  at 
shutter  speed  2  to  1/1000  sec  and  either  TRI-X  black  and  white  print  film  or 
high  speed  Ektachrome  tungsten  color  slide  film. 


2.7.  Operational  Procedures.  Accuracies  and  Flow  Parameters 

The  pressure  and  velocity  measurements  procedure  was  as  follows.  Prior  to 
each  run  with  the  tunnel  running  at  idle,  the  nullity  of  the  pressure  readings 
was  checked.  The  desired  speed  of  the  tunnel  was  set.  The  tunnel  speed, 
extracted  from  the  dynamic  pressure  measurement,  was  checked  against  the 
one  measured  by  laser  Doppler  velocimeter.  The  cavity  gap  was  set  to  the 
desired  width  and  the  frequency  of  cavity  flow  oscillation  was  checked.  Depend¬ 
ing  on  the  degree  of  the  flow  unsteadiness,  the  pressure  averaging  time  was  set 
from  30  to  120  seconds.  The  alignment  of  the  model  was  checked  by  pressure 
reading  differences  of  taps  on  the  four  sides  of  the  downstream  face.  A  cross 
check  of  the  model  alignment  was  made  through  flow  visualization. 

At  the  end  of  a  test  run,  with  the  tunnel  again  set  at  idle,  the  tap  pressures 
were  checked  to  insure  they  had  returned  to  zero.  For  velocity  field  measure¬ 
ments,  the  measuring  volume  was  fixed  at  a  reference  point  on  the  model.  This 
was  done  by  observing  the  position  of  measuring  volume  with  respect  to  the 
model  through  the  eyepiece  on  the  receiving  side  of  the  LDV.  Measurements 
were  made  at  speeds  covering  the  range  of  16  to  27  cm/sec.  The  corresponding 
Re3rnolds  number.  Re,  based  on  the  model  diameter,  d  =  10.16  cm  (4  inches),and 
at  U,  =  23cm Aec  was  24  x  10“^  .  Most  of  the  pressure  data  points  were  meas¬ 
ured  at  least  twice  and  at  some  points  up  to  5  times.  The  repeatability  of  the 
calculated  drag  coefficient  was  about  4%  of  the  average  value. 

2.B.  Reference  Quantities 

Throughout  the  experiment,  the  origin  of  the  reference  coordinate  system 
was  located  at  the  upstream  comer.  The  x  and  y  axes  lie  in  the  symmetry 
plane.  The  x-axis  was  in  the  direction  of  flow;  the  y-axis  was  measured  positive 
outward  with  respect  to  the  center  of  the  cavity.  The  components  of  velocity 
vector  in  the  direction  of  the  x  and  y  axes  were  denoted  u  and  v. 


The  cavity  width,  b,  was  non-dimensionedized  by  normalizing  it  with  respect 
to  Oq  ,  the  momentum  thickness  of  the  upstream  boundary  layer  at  the  point  of 
separation.  Momentum  thickness  for  the  shear  layer  is  defined  as 
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where  £/,  is  the  free  stream  velocity  on  the  outside  edge  of  the  boundary  layer 
at  the  point  of  separation.  Using  £4  instead  of  £/„,  the  free  stream  velocity  in 
front  of  the  model  nose  eliminated  the  need  for  the  blockage  correction.  This 
also  made  the  results  independent  of  the  front  nose  shape  (provided  there 
existed  no  substantial  pressure  gradient  over  the  length  of  the  cavity). 

The  boundary  layer  upstream  of  the  cavity  was  considered  laminar  at 
Ug  =  23cm/sec  and  had  a  shape  factor  of  H  =  — =  2.5  where  6*  =  .1cm.  (dis- 
placement  thickness),  ^  =  0.041cm  and  Ree^^  =  95.  For  the  shear  layer,  the 
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treinsverse  distance  y  was  non-dimensionalized  as  - .  where  is  the 

0 

ZL 

location  of  the  non-dimensional  velocity  -r— with  a  value  of  .5,  and  e  is  the 

Ug 

local  momentum  thickness  of  the  shear  layer  as  defined  by  equation  (2.1).  The 
radius  of  the  axisjmimetric  cavity,  which  is  5  OB  cm  (2  inches)  was  used  to  nor¬ 
malize  the  radial  distance  of  pressure  taps  (r)  on  the  inner  surfaces  of  the  cav- 

ity.  All  drag  forces  were  normalized  by  ,  where  p  is  water  density  and 


A  is  the  cross-sectional  area  of  the  model. 


Chapter  3 


now  Documentation 

The  measurements  of  various  parameters  relevant  to  the  cavity  oscillations 
reported  in  this  chapter  are  not  all  novel,  some  of  them  having  been  reported  in 
references  cited  above  and  in  others.  The  purpose  here  was  to  provide  a  solid 
reference  condition  for  the  measurements  to  be  described  in  Chapters  4  and  5. 
Nevertheless,  some  new  definitions,  clarifications  and  features  are  intro¬ 
duced  . 

3. 1.  Cavity  Flow  Oscillations  as  a  Fiinction  of  b  /e^ 

It  is  known  [Karamcheti  (1956),  Sarohia  (1975),  Rockwell  and  Naudascher 
(1979)]  that  shear  layers  over  cavities  undergo  transition  from  a  steady  layer  to 
an  unsteady  layer  at  a  critical  dimensionless  width  b/Sblmin*  For 
b/00  <  b/®olmiiii  the  shear  layer  bridges  smoothly  over  the  cavity  with  no  distinct 
oscillation  in  it.  The  onset  of  unsteadiness  is  always  associated  with  the  presence 
of  a  distinct  fundamental  frequency,  which  is  amplified  by  the  shear  layer.  In 
this  study,  the  fundamental  frequency  is  identified  as  the  one  that  has  the 
highest  peak  in  the  power  spectrum  of  the  velocity  fluctuations.  Velocity 
fluctuations  were  measured  at  x  =  b  and  y  =1/ .95,(1/ .95  being  the  locus  of 

■jT~=  '  5)  throughout  the  frequency  measurements. 

c/g 

Figure  7  presents  the  variation  of  fundamental  frequency.  F.  with  6/9b  at  a 
fixed  edge  velocity  (f/,  =  33c7n,Aec).  Oscillation  starts  at  6/©o«  81,  which 
represents  the  minimum  cavity  width  6 /©o  I  min.  for  the  present  experimental 
condition.  The  fundamental  frequency  of  oscillation  decreases  as  6/flb 
increases.  This  course  of  frequency  variation  is  interrupted  by  a  sudden  jump 
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of  frequency  from  F  -  4.6  to  6.2  Hz.  This  new  frequency  initiates  another 
course  of  decreasing  frequency  as  b /9b  is  further  increased.  Anon-dimensional 

frequency  can  be  defined  for  these  two  modes  of  oscillation  as  -rr-  (Strouhal 

i/b 

Number).  This  Strouhal  Number  has  a  value  of  «  1  over  the  entire  range  of  the 
steu’ting  mode,  then,  jumps  to  a  value  of  «  1.5  which  it  maintains  over  the  next 
mode.  According  to  Sarohia  (1975)  and  Rockwell  and  Knisely  (1980)  these  two 
modes  correspond  to  Modes  II  and  III  of  the  cavity  fiow  oscillation. 

A  sequence  of  power  spectra  of  the  velocity  fiuctuations  as  a  function  of  b/flb 
is  shown  in  figure  8.  Several  features  of  these  spectra  are  noteworthy.  First  is 
the  presence  of  peaks  both  at  4.6  Hz  and  6.2  Hz  during  the  transition  region 
between  the  modes.  From  these  average  power  spectra  it  is  not  possible  to 
determine  whether  the  two  frequencies  alternate  intermittently  or  whether  both 
coexist  simultaneously.  To  clarify  this  matter,  figure  9  shows  three  oscilloscope 
traces  of  the  velocity  fiuctuations  in  Mode  II,  in  the  transitional  region,  and  in 
Mode  III  of  the  oscillations,  respectively. 

Figure  9b  confirms  the  coexistence  of  two  frequencies  by  showing  a  strong 
modulated  trace  of  the  signal.  A  region  of  hysteresis,  as  shown  in  figure  7,  is 
associated  with  the  transition  region.  Presence  of  frequencies  lower  than  funda¬ 
mental  is  a  prominent  feature  of  Mode  III  (figure  6c).  Such  low  frequencies  have 
been  observed  by  Knisely  (1980)  and  were  shown  to  be  associated  with  the 
interaction  of  successive  shear  layer  vortices  with  the  downstream  corner.  A 
possible  explanation  for  the  existence  of  such  a  low  frequency  is  given  in  the 
flow  model  section  (Section  4.10).  End  of  Mode  III  is  marked  by  a  loss  of  coher¬ 
ence  in  the  oscillation  and  the  disappearance  of  the  fundamental  frequency.  At 
this  stage  {b/6>  160)  .  cavity  flow  begins  to  exhibit  bluff  body  wake  behavior, 
and  a  distinct  frequency  for  this  regime  cannot  be  detected. 
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3.2.  Phase  Heasoranents 

Information  about  the  overall  phase  difference  between  fluctuations  at  the 
upstream  corner  (x  =  0)  and  the  downstream  comer  (x  =  b)  is  needed  for 
improving  our  understanding  of  the  cavity  flow  oscillation.  This  is  the  subject  of 
the  following  discussion. 

The  existence  of  a  strong  phase  gradient  across  laminar  and  turbulent  shear 
layers  is  well  known  [Fre5miuth  (1966),  Sarohia  (1975),  Elder  (1978),  Ziada  and 
Rockwell  (1981),  and  Hussain  (1981)].  This  gradient  can  be  well  approximated  by 
inviscid  theory  (Michalke,  1965),  except  for  regions  close  to  the  downstream 
corner.  Existence  of  such  a  phase  gradient  introduces  difficulties  in  choosing  a 
proper  criterion  to  define  an  overall  phase  difference  between  the  two  corners  of 
the  cavity.  Knisely  (1980),  Ziada  and  Rockwell  (1981),  Rockwell  and  Schachen- 
mann  (1982),  in  their  effort  to  define  an  overall  phase  difference  for  different 
t3T>es  of  self-sustaining  flows,  concluded  that  tracking  values  of  phase  along  the 
edge  of  the  shear  layer  Xu/U,  =  .95)  will  result  in  an  overall  phase  difference  of 
~  2N  TT  .  In  contrast,  other  criteria  ,e.g.,  following  one  of  the  disturbance  ampli¬ 
tude  maxima  within  the  shear  layer,  yield  different  overall  phase  differences  for 
different  geometries  of  self-sustaining  flows. 

In  the  present  study,  phase  measurements  were  obtained  by  following  the 
outer  edge  of  the  shear  layer  along  the  locus  of  u/U,  =  .95.  The  reference  probe 
was  located  at  x  =  b,  y  =  y.gg,  while  the  second  probe  was  moved  along  y^gg  of 
the  shear  layer.  Before  each  set  of  measurements  was  to  be  taken,  the  auto¬ 
correlation  of  the  reference  signal  with  respect  to  itself  or  to  that  of  the  second 
probe  (cross-correlation)  was  obtained.  In  this  way,  any  initial  phase  difference 
due  to  circumferential  or  longitudinal  distances  between  two  probes  was 
detected.  A  sample  of  such  a  measurement  is  given  in  figure  10.  For  the  case  of 
b/<9o  =  85  and  Ug  =  Z3cm/sec  a  typical  phase  angle  distribution  of  the  velocity 


fluctuations  as  a  function  of  x/eo  is  presented  in  figure  11.  The  cavity  shear 

Fb 

layer  has  a  fundamental  frequency  of  5.9  Hz  with  corresponding  77—=  1»  which 

identifies  it  as  the  second  mode  of  oscillations.  Figure  11  depicts  an  overall 
phase  difference  of  47r  between  two  corners  of  the  cavity.  Similar  measurements 
to  determine  the  overall  phase  difference  between  the  two  corners  were  carried 
out  for  the  entire  range  of  Modes  II  and  III  of  cavity  flow  oscillation.  Figure  12 
represents  the  distribution  of  the  overall  phase  difference  as  a  function  of 
It  is  important  to  note  that  the  cavity  flow  retains  a  47t  phase  difference  during 
Mode  11,  while  having  a  Grr  phase  difference  during  Mode  Ill.  This  fact  suggests 
that  overall  phase  difference  between  the  two  corners  be  formulated  as 


2rr 


(3.1) 


where  N  represents  the  mode  of  the  cavity  flow  oscillation.  Using  the  relation 
Uc  =  FX  ,  where  C4  is  the  convection  velocity  (phase  velocity)  of  the  fundamen¬ 
tal  frequency  of  oscillation  and  \  is  the  wavelength,  equation  3.1  can  be  rear¬ 
ranged  as 


2n  X  Uc 
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Equation  3.2  indicates  that  cavity  during  Mode  II  contains  two  wavelengths  of 
the  fundamental  frequency,  while  Mode  III  contains  three  wavelengths.  Also, 
emergence  of  the  Strouhal  Number.based  on  the  phase  velocity  in  equation  3.2, 

yT? 

directs  us  to  define  the  mode  of  the  oscillations  based  on  --7—=  N  rather  than 
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This  convention  defines  each  mode  by  the  number  of  wavelengths 


that  it  contains  and  establishes  a  relation  between  the  mode’s  name  and  the 


corresponding  Strouhal  Number. 

It  is  important  to  note  that  close  to  the  transition  region,  an  important  phase 

deviation  from  2  occurs  (figure  12).  Such  a  deviation  of  overedl  phase 
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difference  edso  was  observed  by  Rockwell  and  Schachenmann  (1982)  for  imping¬ 
ing  turbulent  jets  at  low  Mach  numbers, 

3.3.  Overall  Features  of  Cavity  Flow  Regimes 

3.3.1.  Non-oscillating  Cavity  Flow.  The  non-oscillating  regime  of  the  cavity 
flow  is  marked  by  a  fairly  steady  and  weak  vortex  stationed  inside  of  the  cavity. 
Unsteady  mass  transfer  in  and  out  of  the  cavity  is  minimEd  due  to  lack  of  any 
substantial  oscillations  in  the  shear  layer.  Figure  13  presents  flow  visualization 
of  the  cavity  flow  in  this  regime.  Often,  color  dye  trapped  in  the  cavity  remained 
there  for  a  long  time.  The  effect  of  the  cavity  presence  on  the  outer  flow  field  is 
minimal.  Branching  of  the  shear  layer  at  the  downstream  corner  is  another 
feature  of  this  mode  observed  in  flow  visualization  studies.  It  should  be  men¬ 
tioned  that  for  6/©o  <  6/6^  Imin.^iot-flim  signals  indicated  the  presence  of  weak 
periodic  fluctuations  in  the  shear  layer.  Due  to  lack  of  a  wide  enough  cavity, 
these  perturbations  do  not  amplify  enough  to  make  the  flow  self-sustained.  This 
matter  will  be  discussed  in  more  detail  in  Chapter  5. 

3.3.2.  Self-SvLstaining  Cavity  Flow.  The  definition  of  self-sustaining  flow  and 
the  description  of  the  chain  of  events  that  result  in  flow  self-sustainment  is 
given  in  the  introduction.  As  mentioned  before,  natural  cavity  flow  osciilation  is 
by  definition  self-sustaining. 


3.3.2. 1  General  Characteristics  of  the  Self-sustaining  Modes..  A  self- 
sustaining  oscillation  appears  in  the  cavity  when  6 /So  reaches  6/^  I  min-  Th® 
wavelength  of  the  fundamental  frequency  scales  itself  with  the  width  of  the  cav¬ 
ity  (figure  12).  Flow  in  this  regime  is  stable.  Even  severe  attempts  to  perturb 
the  course  of  oscillations  is  resisted,  depending  on  the  strength  of  the  self- 
sustainment.  In  this  study,  the  most  stable  case  is  defined  as  the  one  that  gen¬ 
erates  the  highest  level  of  organized  fluctuations  for  a  given  b/tfo. 

Overall  phase  difference  between  the  two  corners  of  the  cavity  satisfies  the 

relation  r^=  — =  ^~=N,  where  N  is  the  number  of  wavelengths  contained 
2n  \  Uc 

by  the  cavity  width  during  Nth  mode  of  oscillations.  Shown  in  figure  14  are  flow 
visualizations  of  Modes  II  and  HI  of  the  cavity  flow  oscillations.  In  the  present 
study,  the  cavity  always  started  the  oscillation  with  Mode  II,  omitting  the  first 

mode  y— a  1,  Sarohia  (1975)  observed  similar  omission  of  the  first  mode  in 
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his  study  of  cavity  flow  for  some  free  stream  velocities. 

Instantaneous  flow  visualization  pictures  (figure  14)  show  that  the  wave 
amplitude  grows  as  it  approaches  the  downstream  comer.  This  growth,  which  is 
exponential  according  to  the  linear  theory,  results  in  a  vortex  roll  up.  This  roll 
up  is  very  similar  to  that  observed  by  Brown  and  Roshko  (1974),  Winant  and 
Browand  (1974)  in  the  free  shear  layer.  Mode  HI  shows,  usually,  two  complete 
roll  ups  of  the  shear  layer.  Unlike  the  case  in  a  free  shear  layer,  no  vortex  pair¬ 
ing  is  observed  during  the  cavity  self-sustaining  modes.  Thus,  presence  of  low 
frequency  peaks  in  the  power  spectra  of  Mode  III  of  the  present  study  cannot  be 
associated  with  vortex  coalescence. 

Presence  of  a  strong  stable  vortex  is  another  characteristic  of  the  self- 


sustaining  mode. 
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3.3.2.2  Similarity  to  Externally  Forced  Free  Shear  Layers..  The  ladder-like 
vauiiation  of  cavity  shear  layer  frequency  shows  striking  resemblance  to  the 
externally  forced  tree  shear  layers  (Miksad  1973,  Ho  and  Huang  1982).  One  can 
imagine  that  the  downstream  corner  plays  the  role  of  external  forcing  source. 
Also  the  term  b/[7g  can  be  assumed  to  be  equivalent  to  the  forcing  frequency  in 
an  externally  forced  free  shear  layer  experiment.  Figure  15  shows  the  varia¬ 
tions  of  the  fundamental  cavity  shear  layer  frequency  as  a  function  of  bAJf  In 
this  figure,  several  cases  of  the  cavity  flow  geometries  and  velocities  fall  on  the 
corresponding  mode  lines  of  externally  forced  experiments.  In  general,  one 
must  treat  this  resemblance  with  caution,  because  at  high  subsonic  or  super¬ 
sonic  cavity  flows,  the  wavelength  will  not  scale  with  cavity  width  b.  As  sug¬ 
gested  by  Laufer  and  Monkowltz  (1980)  such  a  scaling  should  follow  the  more 
complicated  feedback  equation 


(3.3) 


where  C  is  the  sound  velocity.  In  the  subsonic  case  C»  U^,  equation  (3.3) 
reduces  to  our  familiar  relation  for  the  Strouhal  Number  7^=  N . 

3.3.3.  Wake  Mode  of  Cavity  Flaw.  The  last  Mode  (III)  of  oscillations,  observed 
in  this  study,  shows  strong  intermittency  as  the  cavity  width  is  increased.  It 
seems  that  eventually  an  intermittency  overcomes  the  coherent  oscillation  and 
the  feedback  mechanism  becomes  ineffective.  The  flow  visualization  studies 
revealed  that  when  the  self-sustaining  regime  ends,  the  shear  layer  no  longer 
reattaches  onto  the  downstream  corner.  The  cavity  flow  becomes  semi- 
periodicaily  unstable  in  a  larger  scale,  comparable  to  the  scale  of  the  model 
(figure  16a,  b).  Time-averaged  flow  visualization  (figure  16c)  shows  that  the 
mean  flow  reattaches  below  the  downstream  corner.  Instantaneous  flow 


■^r’jT'Tirwy 


r 


y: 


26- 


visualization  indicates  that  the  flow  is  highly  asymmetrical  (figure  16a).  The 
flow  field  in  this  regime  shows  some  resemblance  to  a  three-dimensional  blufl 
body  wake  (Werle,1974),  even  though  the  presence  of  the  rear  body  is  still 
strongly  being  felt  by  the  wake  of  the  front  body.  A  detailed  exploration  of  the 
flow  for  cavity-diameters  ratio  fa//?  >  1.65  was  not  carried  out. 

3.4.  Mean  Velocity  Distributions  in  the  Shear  Layer 

For  three  values  of  cavity  width  fa /So  =  85  .fa/flo  =  103  and  b/do  =  130.  mean 
velocity  profiles  were  obtained  at  a  number  of  streamwise  stations.  Edge  velo¬ 
city  (U,)  was  fixed  at  23  cm/sec.  The  cavity  width  b/So  =  65  marks  the  begin- 
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ning  region  of  Mode  II 


■=  2 


.  while  the  cavity  width  fa/(9o  =  103  marks  the 


middle  of  the  same  mode.  For  b/e^  =  130  [77”=  3 
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.  the  cavity  was  oscillating  in 


Mode  III.  Figure  17  presents  nondimensionalized  mean  velocity  profiles  at 
selected  downstream  stations  for  fa/<9o  =  85.  As  shown  in  figure  17,  in  the  early 
stages  of  the  shear  layer  growth,  the  velocity  profile  changes  from  a  boundary 

u 

layer  profile  to  a  shear  layer  profile.  Free-shear-layer  similarity  in  the  77— 

profiles  is  attained  for  fa/9b^  30.  In  figures  18  and  19,  which  present  non¬ 
dimensionalized  mean  velocity  profiles  for  b/Oo  =  103  and  b/Oo  =  130,  respec¬ 
tively,  self-similarity  in  the  -^profiles  is  also  observed 

Figures  20  through  22  present  the  nondimensionalized  mean  transverse  com¬ 
ponent  of  the  velocity,  in  the  shear  layer.  These  figures  reveal  that  unlike 

l4 

the  -^profiles,  mean  transverse  velocity  profiles  do  not  scale  with  the  down- 

f/g 

stream  distance  Thus.no  self-similarity  can  be  found  for  -^profiles.  Gorree, 
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Dimotakis  and  Koochesfahani  (1980)  were  able  to  show  that  for  free  shear 


layers  with  very  high  Re3niolds  numbers, self-similarity  exists  for  -rr— profiles. 

t/g 

The  lack  of  similarity  in  the  -^profiles  of  the  cavity  shear  layers  might  be  due 

L/g 

to  existance  of  the  strong  recirculating  flow  inside  of  the  cavity.  Figure  23 

X 

presents  the  locus  of  values  as  a  function  of  ^  As  in  previous  studies  by 

Browand  (1966),  Sarohia  (1975),  and  Knisely  (1980),  the  y.g  values  were  found 
imtially  to  lie  above  the  level  of  the  separation  edge.  Depending  on  the  stage  of 
the  cavity  flow,  the  locus  of  y.s  approaches  the  y  =0  line  as  the  profile 
develops.  The  shear  layer  shows  a  stronger  deflection  into  the  cavity  for  the 
b/9b  =  103  than  for  b/e  =  85. 

3.5.  Growth  Rate  ol  the  Shear  I^yer 

For  four  values  of  cavity  width,  b/(%  =  85  (Mode  11), 6/^0=  103  (Mode  II). 
b/Oo  =  119  (Mode  II)  and  b/e^  =  135  (Mode  III),  the  momentum  growth  of  the 
shear  layer  is  given  in  figure  24  As  in  the  studies  by  Sarohia  (1975)  and  Knisely 
(1980),  the  momentum  thickness  is  found  to  grow  linearly  with  x/So  Because  of 
this  linear  growth,  one  might  assume  a  self-similarity  for  the  cavity  shear  layers, 
but  it  is  interesting  to  mention  that  in  both  the  previous  studies  and  the 
present  experiment,  the  value  of  the  shear  stress  is  not  constant  (figure  24)  as 
is  required  for  the  self-similarity  of  the  free  shear  layer. 

For  all  cases  with  b/e^^  103,  the  growth  rate  is  fairly  constant  *031 

and  does  not  experience  a  sudden  change  during  the  mode  switching  region.  It 
should  be  noted  that  the  cavity  shear  layer  growth  rate  is  comparable  to  that  of 

turbulent  mixing  layers.  For  -pr— =0,  the  high  Reynolds  number  shear  layer  of 

Liepmann  and  lAufer  ( 1947)  had  a  growth  rate  TtiomenLum  thickness  of  0.035. 
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3.6.  Shear  Stress  and  Velocity  Fluctuation  Distribution  in  the  Shear  Layer 

Measurements  of  shear  stress  and  fluctuating  components  of  the  velocity 
field  are  essential  to  an  understanding  of  the  overall  nature  of  the  cavity  shear 
layer.  As  discussed  previously,  due  to  the  nature  of  the  present  LDV  measure¬ 
ments,  shear  stress  measurements  show  fairly  low  background  noise  levels. 

Such  an  aovantage  was  not  available  for  jr—or  -r— measurements,  (figure  25)  . 
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Figures  26  through  28  depict  the  strearawise  evolution  of  the  transverse 
profile  of  the  shear  stress  for  i)/6b  =  85  (Mode  II).  b/fife  =  103  (Mode  II),  and 

b/0=  130  (III  Mode).  For  large  — ;  shear  stress  profiles  tend  to  show  a  second 

peak  in  the  lower  portion  of  the  shear  layer.  Existence  of  such  a  second  peak 
has  been  observed  [Freymuth(1966) ,Knisely (1980) ] ;  Ziada  and  Rockwell 
(1981)  showed  that  it  can  be  predicted  by  Stuart's  (1967)  theoretical  vortex 
model.  Also  importeuit  for  our  consideration  is  the  fact  that,  even  though  the 

mean  streamwise  velocity  profiles  scale  with  the  maxima  of  the  shear  stress 

00 

profiles  do  not  have  a  constant  value  as  would  be  required  for  self-similarity, 
they  show  an  exponential  growth  with  x/e^  (Figures  29  and  30).  Streamwise 
growth  of  the  maximum  shear  stress  for  several  cavity  widths  is  shovm  in  figures 
29  and  30.  Short  cavity  cases,  b/et  =  85  (Mode  II)  and  b/Bfe  =  103  (Mode  II), 
show  a  region  of  exponential  growth  of  the  maximum  shear  stress.  For  long  cav¬ 
ities  b/6o  >  103,  the  region  of  exponential  growth  is  followed  by  a  saturation 
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region.  In  Mode  II,  the  rate  of  growth  tor  — g — |n^  increases  as  b/6b  increases 

U  g 

(figure  29).  The  transition  region  does  not  show  a  noticeable  change  in  the 
growth  rate,  but  the  general  level  of  the  shear  stress  decreases  for  the  transi¬ 
tion  region  (figure  29). 


•-.  'J 

•  •  x  .  "  .  1 


■  r-.M 

IT.  I 


-29- 


Figure  30  gives  the  distribution  of  maximum  shear  stress  groYrth  for  Mode  III. 
This  mode,  too,  shows  exponential  growth  for  the  maximum  shear  stress  in  the 
shear  layer.  Again.growth  rate  increases  as  6/i9o  increases.  It  is  important  to 
mention  that  the  general  level  of  the  maximum  shear  stress  distribution  is 
higher  for  ft/Sb  =  151.8  and  fe/Sb  =  158.5, which  are  close  to  the  end  of  Mode  III. 

The  maximum  shear  stress  achievable  during  Mode  II,  0.0086,  is  lower  than 
that  in  Mode  III.  which  is  .013.  For  comparison  the  value  of  maximum  shear 
stress  for  a  self-similar,  two-dimensional  turbulent  mixing  layer  as  determined 
by  integrating  the  velocity  profiles  presented  in  Liepmann  and  Laufer  (1947)  is 
.0125.  Integration  of  the  velocity  profiles  measured  by  Wygnanski  and  Fiedler 
(1970)  gives  =  ,015.  Kistler  and  Tan  (1967)  reported  a  maximum 

value  of  0.012  for  two-dimensional  cavity  shear  layers. 

Figures  31  and  32  depict  the  variation  of  the  maximum  streamwise  and 
transverse  fluctuation  as  a  function  of  6/Sb-  Again  an  exponential  growth  of  the 
disturbances  is  evident.  It  is  clear  from  these  figures  that  rms  fluctuations 
suffer  from  a  high  percentage  noise  level  of  either  flow  background  or  electron¬ 
ics. 


3.7.  Mean  and  Iffuctuating  Velocity  Fidd  of  the  Cavity  Flow 

Figures  33  through  37  present  traverses  of  the  velocity  field  through  the  cav- 
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ity  for  one  selected  cavity  width  (6/9b  =  103)  in  Mode  II 
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.  The  particu¬ 


lar  case  selected  has  a  fundamental  frequency  of  5,1  Hz.  Figure  33  presents  the 


mean  streamwise 


u 


velocity  profiles  of  the  upstream  boundary  layer  at  the 


upstream  corner,  at  four  stations  within  the  cavity  and  two  stations  down¬ 


stream  of  the  downstream  corner.  Similarly  figure  34  displays  the  mean 


profiles. 
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The  change  in  flow  direction  inside  the  cavity  in  figures  33  and  34  is  due 
to  the  recirculating  vortex  flow  visible  in  the  flow  pictures.  The  upstream 
boundary  layer  is  laminar  and  has  a  shape  factor  of  2.5.  The  downstream  boun¬ 
dary  layer  profiles  at  x/Bo  =  110  and  x/Oq  =  142  give  respective  boundary  layer 
shape  factors  of  1.54  and  1.46.  These  shape  factors  are  close  to  that  of  a  tur¬ 
bulent  boundary  layer. 

Figures  35  and  36  reveal  that  significant  fluctuations  occur  inside  the  cavity, 
but  the  flow  field  inside  the  cavity  has  a  much  lower  shear  stress  level  than  that 
of  the  shear  layer  (figure  37),  possibly  revealing  that  most  of  the  u'  and  v' 
fluctuations  are  not  well  correlated.  Perhaps  the  most  important  difference 

between  and  jj—  behavior  is  shown  in  the  downstream  boundary  layer 
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profiles  (figure  36).  The  maximum  value  of  the  -j— profiles  in  the  shear  layer 

continues  to  increase,  even  after  the  shear  layer  transforms  to  a  boundary 

layer.  On  the  contrary.the  maximum  value  of  the  profiles  shows  a  decrease 

after  the  downstream  corner.  The  shear  stress  profile  shows  a  monotonic 
decrease  in  the  maximum  value  of  the  profile  as  the  shear  layer  transforms  to 
the  boundary  layer. 


Chapter  4 


The  Drag  Coefficient  of  the  Cavity 

Throughout  this  study  the  term  "drag  coefficient  of  the  cavity"  refers  to  the 
net  force  in  the  flow  direction  being  experienced  by  the  gap.  The  overall  chemges 
in  the  net  drag  coefficient  of  the  modeUwhich  might  arise  from  cavity  presence, 
will  not  be  discussed  in  this  study. 

4.1.  Homentum  Balance  of  the  Cavity 

A  momentum  balance  for  the  fluid  enclosed  by  the  control  surface  shown  in 
sketch  2  gives 

/iZ.(XL- nJtlAi  = -/P  n££A8.4  +  TriLfiAij  (4.1) 

s  s  s 

where  dAi  =  nddx i  ,  dAz,4  =  2nrdr ,  dA^-ndxlx^  and  Jl,  =  ‘ui_+v^  The  x- 


SKETCH  2 


component  of  momentum  balance  will  be 


+  fmd  f  uv  dx  I  =  -Zn  f  |P2  “  P^rdr  +  nd  f  r^dx^  +  iTd,f  Ty^dxg  (4.2) 

Introducing  u=u  +  u' ,v  =v  +«'  to  eq.  (4.2),  taking  the  time  average  and 
neglecting  the  viscous  shear  stress  on  the  free  surface  will  give 

+prr  d  f  tuv  +  irv')dxi  -  “Stt  f  |Pg  -  p’4|rdr  +  ird  f  Ty^  dx^  (4.3) 
Using  the  following  definition  will  reduce  equation  (4.3)  to  equation  (4.4) 


Cdp  =  cavity  pressure  drag  coefficient  =  Zn f^Pz  ~  Pa,  rdr /^pUtird^/^ 


Cos  =  coefficient  of  the  sting  =  -pird^fry^dx^/^U^  rr(iV4 


—pixdjtfiv  +  inr)dxi  =  +  Cop  +  Cos 


(4.4) 


However,  by  definition, the  drag  coefficient  of  the  cavity  can  be  written  as 

Co  -  Cop  +  Cos  (4.5) 

The  combination  of  equations  (4.4)  and  (4.5)  results  in  a  new  relation  for  the 
drag  coefficient  of  the  cavity: 


C])  =  —-nd  f  {mj  +  izV)dzi/'^/w®,  7rdV4 

fit 


(4.6) 


Relation  (4.6)  reveals  that  it  is  possible  to  determine  the  cavity  drag  coefficient 
by  evaluating  the  mean  and  turbulent  momentum  transfer  terms  along  a 
straight  line  connecting  two  corners  of  the  cavity.  Throughout  this  section,  this 
line  will  be  called  y  =  0  line.  Roshko  (1955)  found  that  the  skin  friction 
coefficient  along  the  walls  of  two-dimensional  cavities  were  two  orders  of  magni¬ 
tude  less  than  the  pressure  and  pressure  drag  coefficient.  Therefore,  skin  fric¬ 
tion  in  the  cavity  as  a  contribution  to  the  drag  is  entirely  negligible  and 

Cd  »  Cdp  (4.7) 

4.2.  Measurement  of  the  Cavity  Drag  Coefficient 

The  drag  coefficient  of  the  cavity  is  determined  by  measuring  the  integral  of 
the  pressure  distribution  over  the  solid  walls  of  the  cavity.  Figure  39  presents 
the  variation  of  drag  coefficient  with  b/e^  in  the  three  distinct  regimes  of  cavity 
flow  as  defined  in  chapter  3.  This  figure  reveals  two  rather  significant  charac¬ 
teristics  of  the  cavity  drag.  First,  the  average  level  of  the  drag  coefficient  stays 
low  as  long  as  the  flow  remains  in  either  the  non-oscillating  mode  or  the  self- 
sustained  oscillating  mode.  A  second  feature  is  the  existence  of  a  critical  b/9o, 
in  which,  by  exceeding  its  value,  the  drag  coefficient  of  the  cavity  jumps  to  a 
value  which  is  one  order  of  magnitude  higher  than  the  typical  average  drag 
coefficient  during  oscillating  modes.  During  the  course  of  the  wake  mode  the 
drag  coefficient  increases  monotonically  and  reaches  a  value  of  0.39  at 
b/do  =  250.  No  more  measurements  were  made  beyond  this  point  due  to  model 
limitations. 


4.2.1.  Cavity  Drag  Coefficient  in  Non-OseUlaring  Mode.  The  drag  coefficient 
in  non-oscillating  mode  has  a  typical  value  of  1.3  x  10~*  at  6/(9o  =  77.  When  the 
cavity  was  replaced  by  a  solid  surface,  the  drag  coefficient  of  this  surface  based 
on  the  upstream  boundary  layer  shape  factor  was  estimated  to  be  .7  x  10“®.  In 
turn,  this  drag  coefficient  is  equivalent  to  that  of  6/©o  =  98  of  Mode  11  of  the 
oscillations.  However,  based  on  observations  in  Chapter  3.  the  presence  of  the 
cavity  reduces  the  shape  factor  of  the  downstream  boundary  layer  and  conse¬ 
quently  may  increase  the  skin  friction  drag  on  the  rest  of  the  model.  At  this 
point,  it  is  not  clear  that  a  net  drag  reduction  can  be  obtained  for  the  whole 
model. 


4.2.2.  Exponential  Orowth  of  the  Cavity  Drag  Coefficient  During  Oscillating 
Modes.  Distributions  of  the  drag  coefficient  in  Modes  11  and  111,  plotted  on  a  log- 
linear  scale,  are  presented  in  figure  40.  At  the  beginning  of  the  oscillating 
modes  (Mode  II  in  this  study),  the  drag  coefficient  has  a  value  of  0.9  x  10~®  at 
6/9o  =  65.  The  drag  coefficient  then  increases  exponentially  and  reaches  a  peak 
value  of  0.02  at  b/9o  =  108.  Toward  the  end  of  Mode  II,  the  drag  coefficient 
decreases  and  reaches  a  minimum  value  of  0.004  at  b/9o  =  120.  The  transition 
region  between  Mode  II  and  Mode  III  is  marked  by  an  increase  in  the  drag 
coefficient  value;  Cr)p  =  .0072  at  b/0Q  =  128.  Mode  111  exhibits  another  course 
of  exponential  growth  for  the  drag  coefficient  value.  A  final  value  of  0.39  for  the 
drag  coefficient  is  reached  at  b/e^  =  250  before  transition  to  the  wake 
mode  occurs. 

4.3.  Cavity  Pressure  Distribution 

A  better  understanding  of  the  drag  generation  mechanism  can  be  obtained 
from  pressure  distributions  along  the  solid  surfaces  of  the  cavity. 


Figures  41  through  43  present  the  pressure  coefficient,  'Cp 


H  .  along 


a  radius  of  the  upstream  and  downstream  faces  of  the  cavity  for  several 
selected  values  of  b/e^.  With  each  distribution  the  corresponding  value  of 
is  shown. 

4.3.1.  Pressure  Distribution  for  Non-oscillating  Mode  b/flb  <  82  .  Figure  41a 
reveals  that  in  the  non-oscillating  region,  the  cavity  is  slightly  under  pressure, 
i.e.,  Cp  <  0.  No  significant  differences  were  observed  between  the  upstream  and 
the  downstream  face  pressure  distribution, which  corresponds  to  closed  and 
open  s5mibols,  respectively. 

4.3.2.  Pressure  Dtstribution  for  Oscillating  Modes. 

1)  Mode  II  [82  <  12o) 

Figures  41b  through  41  g  show  the  cavity  pressure  coefficient  distribution  for 
Mode  II  of  the  oscillation.  All  cases  have  negative  pressure  distribution  over  the 
entire  upstream  face.  The  most  important  point  to  note  in  all  cases  is  the  steep 
pressure  rise  near  the  outer  edge  of  the  downstream  face.  In  each  case,  Cp 
reaches  its  maximum  positive  value  {Cp^^  at  the  corner.  The  magnitude  of 
Ca  reaches  a  maximum  at  b/Sn  =  108.  then  decreases  toward  the  end  of  Mode 
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11.  The  downstream  face  experiences  a  negative  pressure  below  the  downstream 
comer  with  the  same  level  and  trend  as  that  of  the  upstream  face.  Thus,  the 
contribution  to  the  total  drag  from  these  two  portions  of  the  cavity  walls  cancel 
each  other.  This  leaves  the  outer  edge  of  the  downstream  face  with  >  0  ,  tbe 
main  contributor  to  the  cavity  drag.  It  is  also  interesting  to  mention  that  in  all 
cases,  the  average  pressure  on  the  upstream  face  is  negative;  and  thus  this  face 
experiences  thrust  during  Mode  II  of  the  oscillations. 

The  pressure  buildup  near  the  upstream  corner  is  evidently  due  to  the 
deflection  of  the  shear  layer  into  the  cavity.  This  deflection  of  the  shear  layer 


causes  a  strong  inflow  to  the  cavity  and  contributes  to  the  generation  of  a  stable 
vortex.  A  large  portion  of  the  cavity  experiences  negative  pressure  coefficient 
due  to  the  presence  of  strong  recirculating  flow  inside  of  the  cavity.  Similar 
types  of  observations  were  made  by  Johannesen(1955) ,  Roshko(1955),  McGregor 
and  YThite  (1970),  Charwat,  et  al.  (1961),  Fox  (1964)  and  Rossiter  (1964). 

2)  Mode  III  120^  6/9o^  155 

Figures  42a  through  42d  present  pressure  distribution  for  Mode  III.  Except 
for  a  smoother  pressure  rise  near  the  outer  edge,  the  pressure  distribution  on 
the  upstream  face  shows  behavior  similar  to  that  of  Mode  II. 

4.3.3.  Pressure  Distribution  for  the  Wake  Mode.  Pressure  distributions  for 
three  cases  of  the  wake  mode  are  shown  in  figure  43  .  Two  main  features  dis¬ 
tinguish  the  wake  mode  from  the  oscillating  mode.  First,  the  entire  upstream 
face  of  the  cavity  experiences  positive  pressure,  i.e.,  Cp  >  0.  Second,  the  stagna¬ 
tion  point,  the  position  of  Cp^^,  steadily  migrates  toward  the  center  of  the  cav¬ 
ity.  We  expect  that  for  a  large  enough  value  of  b,  the  stagnation  point  should 
settle  at  the  center  of  downstream  face.  While  pressure  on  the  rear  face  is 
increasing,  the  average  base  pressure  level  decreases  as  b/fib  increases  (figure 

43  ). 

4.3.4.  Overall  Characteristics  of  the  Cavity  Pressure  Distribution.  Figure  44 
presents  a  collection  of  typical  upstream  corner  pressure  distributions  for 
different  regimes  of  the  cavity  flow.  The  start  of  the  self-sustained  oscillations  is 
accompanied  by  an  overall  drop  of  the  upstream  face  pressure  distribution.  The 
transition  from  Mode  II  to  Mode  III  (6/6^  =  126)  is  marked  by  a  substantial 
increase  in  upstream  face  pressure.  A  flow  visualization  study  revealed 
that  the  absence  of  a  '^rong  recirculating  flow  inside  the  cavity  might  be  respon¬ 
sible  for  the  pressure  rise  of  the  upstream  face.  In  Mode  III,  low  upstream  face 
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pressure  recovers  and  shows  a  level  similar  to  that  of  Mode  11. 

Upon  the  loss  of  the  self -sustained  oscillations,  another  dreistic  decrease  in 
the  average  level  of  the  upstream  face  pressure  distribution  occurs.  At 
6/(9b  =  210  [b/R  w  1.7),  the  average  upstream  face  pressure  coefficient  has  a 
value  of  -0.13. 


«  •  •  ^ 


Figure  45  represents  variation  of  maximum  pressure.  Cp^,  with  b/(^.  The 
value  of  Cj,  increases  as  6 /Sb  increases  until  it  reaches  a  local  maximum 
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at  6/0b  =  108.  This  maximum  is  immediately  followed  by  a  valley  which  marks 
the  last  point  of  Mode  II.  The  transition  region  and  Mode  III  show  a  monotonic 
increase  in  the  value  of  Cp^^.  Appearance  of  the  wake  mode  is  marked  by  a  sud¬ 
den  increase  in  the  value  of  CL  from  +0.17  to  +0.43  at  b/9n  w  162.  The  value 
of  shows  slight  decrease  during  the  wake  mode.  The  postion  of  the  stagna¬ 
tion  point  shifts  toweird  the  center  of  the  cavity  as  b/efe  increases.  It  is  impor¬ 
tant  that  the  location  of  is  fixed  in  the  vicinity  of  the  upstream  corner  in 
the  self -sustained  modes.  Migration  of  this  point  toward  the  center  of  the  cavity 
starts  only  when  the  cavity  loses  the  self -sustained  oscillation  to  the  wake  mode. 
The  rather  dramatic  change  in  the  pressure  field  of  the  cavity  caused  by  exceed¬ 
ing  a  critical  value  of  b/e^  has  been  observed  by  Roshko  (1955),  Charwat  (1961), 
Fox  (1968)  and  Maul  and  East  (1963).  However,  the  role  of  self-sustained 
oscillations  in  determining  this  critical  b/0p  was  in  general  unnoticed. 

Variation  of  Cup  at  the  middle  of  two-dimensional  cavities  studied  by  Roshko 
(1955)  and  Fox  (1968)  as  well  as  those  of  the  present  study  are  shown  In 
fig.  46.  It  is  Interesting  that  in  all  three  cases  the  location  of  steep 
jump  in  the  value  of  C  agrees  fairly  well. 
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4.4.  Observation  Sadden  Drag  and  Pressnre  Change  in  Other  Flows 

Several  other  investigators  have  reported  a  sudden  change  from  a  steady  low- 
drag  flow  to  a  highly  non-steady  high-dreig  flow  over  cavities,  cutouts,  or  bluff 
bodies  in  tandem.  The  intention  in  this  section  is  to  compare  some  of  their 
results  with  those  of  the  present  study,  and  to  give  a  possible  explanation  for 
the  existence  of  the  critical  length,  b/&Q, 

Koenig  (1978)  in  his  study  of  flow  over  the  combination  of  a  disk  and  a 
cylinder,  defined  two  regimes  based  on  the  drag  behavior  of  the  ss^stem,  neunely, 
sub-  and  super-critical.  Drag  in  a  subcritical  regime  has  values  one  order  of 
magnitude  lower  than  that  of  a  supercritical  regime.  Separated  shear  layer  in  a 
subcritical  regime  smoothly  joins  onto  the  rear  body,  while  in  the  supercritical 
regime  the  shear  layer  becomes  unsteady  on  a  much  greater  scale.  Koenig's  flow 
visualization  study  revealed  that  the  optimum  low  drag  case  appears  to  be  a 
more  general  example  of  flow  over  cavities.  Figure  47  presents  optimum  drag 
measurements  (Koenig,  1978)  superimposed  on  the  present  drag  measurements. 
Despite  the  basic  differences  in  the  geometry  and  the  curvature  of  the  shear 
layer,  figure  47  confirms  Koenig’s  flow  visualization  conclusions. 

Koenig  ailso  reports  the  presence  of  a  steady  and  stable  flow  field  in  the 
supercritical  region.  Production  of  sharp  acoustic  tones  was  another  feature  of 
subcritical  regime  that  was  reported  by  Koenig.  The  symptoms  of  subcritical 
regime  are  similar  to  those  of  simple»self-sustained.cavity  flow  oscillations. 

Mair  (1965)  reported  a  sudden  change  from  steady  low-drag  flow  to  nonsteady 
high-drag  flow  near  the  base  of  a  streamwise  oriented  cylinder  as  the  gap 
between  the  cylinder  and  a  downstream  disk  was  changed.  Zdravkovich  (1977) 
observed  an  abrupt  start  of  vortex  shedding  from  an  upstream  cylinder  accom¬ 
panied  by  an  increase  in  drag  of  a  S3^tem  comprised  of  two  cylinders  in  tandem 
when  spacing  between  them  was  changed. 


The  presence  of  a  feedback  loop  actively  doing  the  task  of  cross¬ 
communication  between  two  corners  in  the  gap  generated  by  two  bluff  bodies  in 
tandem  may  possibly  explain  extra-stability  and  low-drag  characteristics  of  such 
flow  systems. 

The  lock-on  of  the  separated  shear  layer,  or  the  position  of  Cp^  to  the 
downstream  corner  due  to  feedback  mechanism,  has  been  overlooked  by  previ¬ 
ous  investigators  and  has  often  been  considered  a  secondary  effect.  A  further 
discussion  of  this  matter  will  be  given  in  section  4.10. 

4.5.  Estimate  of  Co  by  Evaluation  of  ^-and  Terms  Along  y  =  0  line 

An  estimate  of  Cp  can  be  obtained  by  integrating  the  shear  stress  and  mean 
momentum  term  over  the  periphery  of  the  cavity  (equation  4.6).  LDV  technique 
was  used  to  measure  the  x  and  y  component  of  the  velocity  vector,  to  determine 

M_and  terms  along  y  =  0  line.  Tor  each  case  between  20  to  30,  measure- 
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ment  points  were  selected  along  the  y  =  0  line,  depending  on  cavity  length. 
Smaller  spacing  between  measurement  points  was  selected  close  to  the  down¬ 
stream  corner.  Due  to  the  flow  field  symmetry  of  the  wake  mode,  no  repeatable 
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measurements  could  be  done  in  that  regune.  Thus,  evaluation  of  and 
terms  were  limited  to  the  oscillating  regime  of  the  cavity  flow. 

4.5.1.  ConfiTmation  of  Exponentiiil  Orowth  far  Cp.  Figure  40  presents  the 
variation  of  Cp  based  on  the  integration  of  — g— and  along  the  y  -  0  line 

I/q  Ug 

with  ft /9b-  It  is  important  to  note  that  the  exponential  growth  of  Cp  during  oscil¬ 
lating  modes  has  been  confirmed  by  this  independent  method.  Except  for  the 
cases  close  to  the  end  of  Mode  II  and  close  to  the  transition  region  between 
Modes  II  and  111,  the  absolute  values  of  Cp  and  Cp  agree  fairly  well.  Study  of 
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the  77— .  77-:  ^  euid  — =-  distributions  along  the  y  =  0  line  is  essential  in 

t/,  t/,  £/,  U; 

further  understanding  of  the  mechanism  which  is  responsible  for  the  Cp 
behavior.  This  is  the  subject  of  the  following  chapters. 

4.6.  Distribution  of  Along  y  =  0  line 

Provided  that  the  distribution  of  €>  at  a  fixed  x  position  is  uniform  on  the 
periphery  of  the  cavity,  the  integral  of  v  along  the  y  =  0  line  should  be  minimal 
for  a  satisfactory  mass  balance  of  the  cavity. 

Starting  from  figures  48a  through  48g,  a  sequence  of  -j—  distributions  for 

a  few  values  of  &/0o  in  Mode  II  is  given.  Each  figure  is  accompanied  by  its 
corresponding  dreig  coefficient. 

A  common  feature  of  all  77— distributions  is  the  existence  of  a  narrow  region 

C/, 

of  inflow  to  the  cavity  -v  close  to  the  downstream  corner.  The  negative  v 
region  is  balanced  by  a  broader  region  of  the  outflow  from  the  cavity  +17  close  to 
the  upstream  corner.  The  stronger  the  inshooting  of  flow  into  the  cavity,  the 
higher  is  the  peak  of  the  positive  v  region.  The  peak  of  H7  increases  as  h/e^ 
increases  and  reaches  a  maximum  at  b/e-  112.5,  which  corresponds  to  the 
maximum  drag  case.  This  peak  eventually  decreases  toward  the  end  of  Mode  II. 
In  Mode  III  the  i7  region  becomes  slightly  broader  and  shows  a  small  change 
(figure  49). 

V 

In  terms  of  stable  vortex  inside  the  cavity,  — — ;  distribution  can  be  inter- 

preted  as  the  induced  velocity  of  cavity  vortex  on  the  y  =  0  line.  Therefore,  it 
can  be  concluded  that  the  strength  of  the  vortex  increases  during  the  Mode  II  as 
6/9o  increases.  The  strongest  and  most  stable  vortex  appears  at  the  drag 


maximum  case;  i.e.,  at  b/0o  =  11S.5.  F  =  5  Hz.  The  cavity  vortex  becomes  weedcer 
during  transition  and  recovers  its  strength  during  Mode  III. 


A  check  of  mass  balance  showed  that  between  93  and  97  percent  of  mass 
entrance,  exits  the  cavity.  Either  a  slight  angle  of  attack  or 

presence  of  three-dimensionalities  in  flow  are  responsible  for  imperfecLions  in 
the  mass  balance.  It  remains  to  mention  that  it  is  logical  to  assume  that  positive 
entrainment  to  the  shear  layer,  at  the  upstream  region,  is  the  primary  process 
for  generating  +v  region  and  consequently  causing  the  deflection  of  the  shear 
layer  into  the  cavity  at  the  downstream  corner.  Relation  between  the  positive 
entrainment  and  amplification  of  the  perturbations  in  the  initial  region  of  the 
shear  layer  will  be  discussed  in  the  flow  model  section. 


4.7.  Distribution  of  jj—  Along  y  =  0  line 

The  -“distribution  along  y  =  0  line  (figures  50a  through  50g)  exhibits  some 
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interesting  features.  A  t5T5ical  distribution  shows  that  -—increases  from  zero 

^0 

at  X  =  0  to  a  peak  value,  and  is  usually  followed  by  a  stagnation  flow-like 
behavior  where  velocity  experiences  a  sharp  decrease  down  to  zero  at  x  =  b. 

The  stagnation  region  often  shows  several  local  peaks  and  valleys.  In  Mode  II, 


the  peak  value  of 


shows  an  increase  until  it  reaches  its  maximum  at 


6 /9c  =  112.5  ,  corresponding  to  the  maximum  drag  case.  The  end  of  Mode  II  and 
the  transitional  region  are  marked  by  lower  values  for  the  peak  of  similar 

to  that  of  the  distribution. 


Mode  III  does  not  show  a  substantial  change  in 


distribution  or  value  of 


U, 

peak  velocity  during  its  course  (figure  51).  It  is  important  to  note  that  -Tr—imax 


changes  from  0.3  to  0.6  during  its  rise  in  Mode  II.  which  is  a  cleeu"  indication  of 
defiection  of  shear  layer  into  the  cavity.  Consequently  this  deflection  will  cause 
the  stagnation  of  the  high  velocity  portion  of  the  shear  layer  velocity  profile  at 
the  downstream  corner.  It  seems  that  the  natural  diffusion  of  the  shear  layer 
due  to  the  cavity  width  increase  (shear  layer  length),  plays  a  less  important  role 
in  bringing  higher  velocity  fiuid  into  the  stagnation  region.  Decrease  in  the  peak 

value  of  during  the  transitional  region,  despite  the  increase  of  cavity 

L/fl 


width,  supports  the  previous  argument.  In  Mode  111,  the  averaged  level  of  the 

-^Imax  (while  having  a  value  of  0.5)  stays  constant  (figure  51).  Diffusion  of  the 

shear  layer  into  the  cavity  in  Mode  III  probably  plays  a  more  important  role  than 
the  defiection  of  the  shear  layer  in  bringing  high  velocity  fluid  to  the  stagnation 
region. 


4.6.  Evaluation  of 


-7^  Along  y  =  0  Line 


£// 


The  -^j^term  is  eveduated  along  the  y  =  0  line,  using  the  independent  meas¬ 


urements  of  y^and  and  the  results  are  shown  in  figures  52a  through  52g. 


u. 
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A  typical  distribution  of  shows  a  wide  region  of  positive  values  for 

Ub  Ub 


which  is  followed  by  a  sharp  negative  peak  close  to  the  downstream  corner.  The 
region  of  positive  values  of  contributes  thrust  to  the  net  force  balance  while 


the  region  of  negative  values  of  contributes  drag,  (figure  53) . 


is  £ilwa3rs  negative  during  Modes  II  and  III  of  the 


It  was  found  that  f 
0 

oscillations.  Even  though  -^distribution  showed  a  feiirly  good  mass  balance 

for  the  cavity,  the  -C  region  was  weighted  by  higher  values  of  -MZ  during 
integration  along  the  y  =  0  line,  resulting  in  a  positive  net  drag  for  the  cavity. 


For  convenience,  the  following  conventions  are  introduced.  The  net  contribu- 
tion  of  — 5- term  to  total  drag  coefficient 


Net  contribution  of 


— term  to  total  drag  coefficient 

ui 


-  ^Dtnr 


Total  drag  coefficient  =  Cq. 
Therefore 


~  Qtbp  Qwtt 

A  comparison  between  QjOP  and  Qwv  contributions  to  the  net  drag  will  be 
given  in  the  next  section. 


4.9.  Shear  Stress  Distribatioii'Alaas  y=0  line 


Figures  54a  through  54g  present  shear  stress  distribution  for  Mode  II.  Shear 
stress  grows  along  the  y  =  0  line,  reaches  a  peak  and  vanishes  steeply  as  the  flow 
approaches  the  downstream  stagnation  point.  Often  two  peaks  appeeir  before 

the  decrease  in  the  stagnation  region  starts.  The  peak  of  reaches  the  max¬ 


imum  value  of  .0063  at  6 /fife  =  112.5,  which  corresponds  to  the  maximum  drag 
case.  During  the  transition  region  the  peak  value  of  the  shear  stress  decreases 
to  0.003. 


The  start  of  Mode  III  is  marked  by  the  appearance  of  a  wide  region  of  high 
shear  stress  value  close  to  the  downstream  comer.  In  Mode  III  (figure  55),  the 

peak  value  of  increases  from  .007  at  b/9o  =  132  to  0.0098  at  b/Bg  =  146. 
Uf 


Figure  56  presents  the  variation  of 


— —  euid  — - — versus  6 /fife.  It  also  reveals 


that  in  Mode  II,  especially  in  the  beginning  stages,  contributes  up  to  80%  of 
the  total  drag  coefficient.  As  b/Sfe  increases,  the  role  of  becomes  more 

important  but  its  contribution  never  exceeds  that  of  ^Dao-  Mode  111  shows  a 
dominant  role  for  Cj^^.  Immediately  after  transition,  C/jmr  constitutes  60%  of 
Cj)  and  increases  to  70%  at  the  critical  gap  {b/Oo  =  146).  The  message  of  figure 

56  is  that  consideration  of  the  shear  stress,  as  the  only  means  of  momen¬ 
ta. 


turn  transfer  along  the  cavity  y  =  0  line  will  result  in  a  30  to  80%  underestima¬ 
tion  of  the  cavity  drag  depending  on  cavity  length. 


Kistler  and  Tan  (1967)  made  such  estimations  oE  the  drag  coefficient  of  a 
two-dimensional  cavity  with  a  turbulent  upstream  boundary  layer.  Using  hot 
wire  anemometry,  he  determined  the  shear  stress  distribution  along  the  cavity 
mouth.  His  results  show  up  to  60%  underestimation  of  Cd  with  respect  to  Cop. 


It  is  important  to  recall  that,  due  to  insensitivity  of  hot-wire  techniques  to  the 

fZv 

direction  of  the  flow,  it  is  impossible  to  accurately  estimate  the  term  along 
the  y  =  0  line. 

4.10.  A  Possible  now  Model 

The  drag  coefficient  of  the  cavity  stays  low  as  long  as  self-sustained  oscilla¬ 
tions  are  present  in  the  flow.  Flow  visualization  studies  have  revealed  that  the 
shear  layer  (time  averaged)  re-attaches  onto  the  downstream  corner.  Only  dur¬ 
ing  the  wake  mode  does  the  shear  layer  stagnate  below  the  downstream  corner. 
The  fact  that  the  stagnation  point  (average  position  of  maximum  pressure)  is 
always  fixed  at  the  downstream  corner  supports  the  fiow  visualization  conclu¬ 
sion.  There  remains  one  question  to  answer:  what  mechanism  keeps  the  flow 
from  stagnating  below  the  downstream  corner  and  consequently  results  in  a  low 
drag  value  for  the  cavity?  What  follows  is  a  possible  answer  to  the  above  ques¬ 
tion. 

It  is  logical  to  assume  that  the  volume  of  flow  that  is  entrained  by  the  shear 
layer  is  proportional  to  the  strength  of  the  shear  forces  in  the  shear  layer.  Fig¬ 
ure  57  reveals  that  the  amount  of  fluid  that  is  entrained  by  the  shear  layer  (or 

exits  the  cavity  v±c  |+.„j  follows  the  trend  of  maximum  shear  stress  of  the 
shear  layer  at  a  fixed  streeimwise  position;  e.g.,  x/9o  =  65. 

On  the  other  hand,  mass  balance  of  the  cavity  requires  that  the  shear  layer 
deflect  into  the  cavity  in  order  to  supply  fluid  for  positive  entrainment  to  its  ini¬ 
tial  region.  Figure  23  presents  the  location  of  7—=  0.5  in  the  shear  layer.  At 
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the  beginning  of  Mode  II  (6/9o  =  83,  F  =  6  ).  where  shear  stress  level  and  growth 
is  low,  the  shear  layer  displays  mild  curvature.  At  the  maximum  drag  case 
(6/9o  =  1 12.5),  where  the  shear  stress  level  throughout  the  shear  layer  is 


higher,  the  shear  layer  shows  steeper  curvature  at  the  initial  region  and  a  sharp 
deflection  into  the  cavity.  Therefore,  the  degree  of  the  deflection  of  the  shear 
layer  depends  on  the  nature  of  the  instability  growth  in  the  shear  layer.  The 
severity  of  the  shear  layer  deflection  is  decisive  in  the  nature  of  the  shear  layer 
interaction  with  the  downstream  corner. 

Rockwell  and  Knisely  (1979)  show  that  if  the  respective  location  of  the  down¬ 
stream  comer  is  raised  by  -^=  6,  the  rms  transverse  fluctuations,  iT,  will 

So 

drop  by  a  factor  of  2.  Sarohia  (1976)  demonstrated  that  by  external  mass  injec¬ 
tion  to  the  cavity,  it  is  possible  to  deflect  back  the  shear  layer  at  the  corner  and 
weaken  the  strength  of  the  cavity  flow  fluctuations.  Based  on  these  experiments 
it  can  be  concluded  that  for  given  b/e^  the  displacement  of  the  impinging  shear 
layer  on  the  corner  due  to  deflection  should  occur  in  such  a  way  that  the  feed¬ 
back  loop  would  not  lose  its  strength.  Any  severe  external  interference  with  the 
established  self-sustained  oscillations  would  result  in  an  overall  drop  of  the 
fluctuation  level  throughout  the  shear  layer.  These  two  experiments  provide  the 
missing  link  between  the  role  of  the  corner  and  the  fluctuation  level.  Once  the 
feedback  loop  due  to  the  frequency  selection  process  is  established,  any  inboard 
displacement  of  the  velocity  profile  caused  by  the  natural  tendency  of  the  wake 
mode  will  be  resisted  by  the  folloyring  mechanism.  A  negative  displacement  of 
velocity  profile  at  the  downstream  corner  will  result  in  the  lowering  of  the  shear 
stress  level  which  will  be  followed  by  a  decrease  in  the  volume  of  positive 
entrainment.  However,  a  less  positive  entrainment  would  result  in  reduced 
deflection  of  the  shear  layer  and  restoration  of  the  initial  position  of  the  shear 
layer  velocity  profile  at  the  downstream  corner.  This  mechanism  will  not  allow 
the  full  velocity  profile  at  the  impinging  corner  to  be  displaced  more  than  its 
total  thickness,  and  as  a  consequence,  keeps  the  location  of  flxed  at  the 

downstream  corner.  The  strength  of  this  restoration  mechanism  depends  on 
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the  stability  of  the  cavity  flow  oscillations,  which  was  already  established  by  the 
frequency  selection  process.  It  is  possible  that  because  of  the  existence  of  such 
a  restoration  mechanism,  some  low  frequency  oscillations  appear  in  the  shear 
layer.  According  to  the  above  process  these  frequencies  need  not  be  the  har¬ 
monics  of  fundamental  frequency.  As  stated,  Knisely  (1980)  observed 
0.5,  0.4  and  0.6  of  fundamental  frequency  in  the  shear  layer  over  cavities.  In 
the  present  study  severe  low  frequency  components  were  observed  in  Mode  III 
(last  mode)  of  the  oscillations.  Toward  the  end  of  Mode  III,  as  a  result  of  a 
decrease  in  the  effectiveness  of  the  feedback  loop,  the  restoration  force  became 
weaker.  Consequently,  the  natural  wake  mode  of  the  front  body  overcame  the 
restoration  mechanism,  and  the  shear  layer  had  more  freedom  to  move  below 
the  downstream  corner  and  stagnate  inside  the  cavity. 


Chapter  5 


Control  of  Cavity  Roir  OocillationB 

One  of  the  objectives  of  the  present  study  was  to  develop  and  use  a  technique 
which  is  capable  of  controlling  the  frequency  and  amplitude  of  established  flow 
oscillations.  However,  such  a  technique  can  be  used  to  obtain  information 
regarding  the  necessary  conditions  for  the  onset  of  self -sustained  oscillations. 

This  chapter  will  report  on  emd  discuss  a  series  of  experiments  regarding  the 
application  of  a  new  technique  for  forcing  the  cavity  shear  layer.  Throughout 
these  experiments,  the  response  of  the  cavity  shear  layer,  both  in  non¬ 
oscillating  and  oscillating  modes,  is  studied  against  different  forcing  frequencies 
and  amplitudes. 

5.1.  Forcing  of  the  Shear  Layers 

Forcing  of  the  unstable  and  stable  modes  of  laminar  or  turbulent  shear 
layers  has  been  the  topic  of  many  previous  studies.  Basically,  forcing  techniques 
fall  into  two  catagories: 

1)  Perturbations  to  the  bulk  of  either  flow  stream  (Miksad,  1972,  and  Ho  and 
Huang,  1982). 

2)  Localized  perturbations  in  the  initial  region  of  the  shear  layer  (Oster  and 
Wygnanski,  1982). 

In  the  first  t)q)e  of  forcing,  e.g.,  by  oscillating  one  of  the  streams,  the  whole 
streamwise  length  of  the  shear  layer  experiences  forcing  effects.  As  a  result  no 
accurate  information  can  be  deduced  about  the  influence  of  perturbations  on 
the  initial  regions  of  the  shear  layer. 
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Wygnanskl  and  Osier  showed  that  forcing  effects  can  be  localized  by 
sinusoidal  oscillation  of  a  pivoted  thin  flap  at  the  leading  edge  of  the  shear  layer 
splitter  plate.  The  other  technique  used  to  localize  the  effect  of  forcing  is  to 
introduce  controlled  excited  (T-S)*  waves,  amplified  by  the  boundary  layer,  into 
the  shear  layer.  The  idea  of  actively  forcing  the  boundary  layer  was  first  pro¬ 
posed  by  Liepmann  (1943).  The  most  popular  technique  to  facilitate  the  idea 
was  the  vibrating-ribbon  developed  by  Schubauer  and  Skramstad  (1947).  The 
difficulty  with  either  the  vibrating-ribbon  or  the  flapping  edge  technique  for  the 
present  experiments  is  that  it  cannot  be  adapted  to  axis)niimetric  flows. 

5.2.  Surface  Heating  Technique 

Liepmann,  Brown  and  Nosenchuck  (1982)  recently  showed  that  time- 
dependent  surface  heating  cam  be  used  to  introduce  perturbations  into  laminar 
boundary  layers.  Such  a  technique  can  be  easily  adapted  to  excite  T-S  waves  in 
an  axis3mimetric  boundary  layer  such  as  that  of  the  present  study.  The  idea  of 
periodic  heat  release  into  the  boundary  layer  for  the  purpose  of  exciting  T-S 
waves  was  also  realized  by  Kendall  (1967 )  in  supersonic  boundary  layers.  By 
appl3dng  periodic  current  to  two  closely  spaced  parallel  wires  flush  mounted  on 
the  wall,  Kendall  was  able  to  initiate  glowing  between  the  two  wires  and  conse¬ 
quent  heat  release  into  the  boundary  layer.  As  a  result  of  this  periodic  heat 
release  into  the  boundary  layer.  T-S  waves  of  up  to  30  kHz  were  excited. 
Requirement  of  low  pressure  to  attain  glowing  makes  this  technique  of  limited 
use  for  subsonic  flow  applications. 

For  cases  in  water,  Liepmann,  et  al.  (1982)  show  that  the  effect  of  "heating 
€Uid  cooling  of  the  surface  is  roughly  equivalent  to  suction  and  blowing  or  nega¬ 
tive  and  positive  surface  displacements."  with 
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where  fi  is  dynamic  viscosity.  T  is  the  temperature  of  the  fluid,  Pr  is  the 
Prandtl  number,  w  is  the  width  of  the  strip,  and  /?,  is  the  Reynolds  number. 
Such  an  equivalent  normal  surface  oscillation  would  periodically  redistribute 
vorticity.  Depending  on  the  state  of  the  boundary  layer,  the  external  perturba¬ 
tions  would  either  be  dampened  or  amplified  exponentially.  Linear  theory  can 
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predict  the  exponential  growth  of  the  (T-S)  waves,  based  on  the  parameters  -jr— 
and  Red*  for  a  given  pressure  distribution. 

5.2.1.  v^piicafioTi  of  T-S  Waves  as  a  Means  of  Forcing  of  the  Shear  Layer. 
Forcing  T-S  waves  in  boundary  layers  has  been  practiced  for  a  long  time.  How¬ 
ever,  to  our  knowledge,  no  use  has  been  made  of  T-S  waves  to  force  a  separated 
shear  layer  resulting  from  such  a  boundary  layer.  In  the  present  study, 
inspired  by  the  experiment  of  Uepmann,  et  al..  periodic  heating  of  a  thin  film 
was  used  to  excite  the  T-S  waves  in  the  boundary  layer  upstream  of  the  cavity 
shear  layer,  benefltting  from  amplification  of  the  waves  before  their  introduc¬ 
tion  into  the  shear  layer  (figure  58).  One  can  imagine  that  the  effect  of  the  T-S 
waves  or  v"  component  at  the  upstream  corner  of  the  cavity  would  be  similar  to 
that  of  the  oscillating  flap  at  the  leading  edge  of  the  splitter  plate  in  Oster  and 
Wygnanski’s  experiment.  It  will  be  seen  later  that  the  frequencies  needed  to 
excite  the  free  shear  layer  are  available  in  the  range  of  boundary  layer  T-S 
waves.  Thus,  a  transformation  of  boundary  layer  T-S  waves  to  the  shear  layer 
Kelvin-Helmholtz  waves  is  possible. 


5.2.2.  Strip  Heater  Set-up.  Initially  three  strip  heaters  were  mounted  on  the 
ellipsoidal  nose  of  the  model  (figure  59)  at  s  =  9.2  cm.  12  cm,  and  15.5  cm  where 
s  is  the  distance  on  the  model  surface  from  the  nose  of  the  model.  Double  stick 
scotch  tape  was  used  to  electrically  isolate  heaters  from  the  model.  The  first 
and  third  heaters  were  stainless  steel  ribbons  3.7  mm  (.125  inch)  wide,  .051  mm 
(.002  inch)  thick,  wrapped  around  the  front  nose  of  the  model.  The  second 
heater  was  narrower.  ,08  cm  (.2  inch)  wide.  Each  heater  entered  the  model  from 
a  flush  mounted  round  adaptor,  designed  to  connect  the  heaters  to  the  current 
supply  lines.  The  resistance  of  the  strip  heaters  was  2  -  4  ohms. 

A  high  quality  commercial  dual  channel  power  amplifier  (Hafler  DH-500)  was 
used  to  supply  current  to  the  strip  heaters.  DH-500  is  a  500  watts  per  channel 
amplifier  which  can  be  easily  modified  to  provide  frequencies  even  below  1  Hz 
while  delivering  adequate  power  to  the  strip  heater.  An  HP  Function  Generator 
(3325A)  was  used  to  provide  the  selected  frequency  to  the  power  amplifier.  It 
should  be  mentioned  that  due  to  the  quadratic  dependence  of  Joule  heating  on 
the  voltage,  the  input  frequency  to  the  power  amplifier  had  to  be  half  of  the 
desired  forcing  frequency. 

5.2.3.  Important  Parameters  of  Strip  Heater  Technique.  Several  parameters 
are  important  for  the  efllciency  of  a  strip  heater  in  forcing  both  the  T-S  waves 
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and  the  shear  layer.  Uepmann,  et  al.  (1982)  stated  (where  W  is  the  strip 
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heater  width  and  6*  is  the  boundary  layer  displacement  thickness)  is  an  impor¬ 
tant  parameter  in  effectiveness  of  strip  heater  operation  on  the  local  forcing  of 
T-S  waves.  The  state  of  the  free  stream  pressure  gradient  is  another  important 

W 

factor.  In  the  present  study,  the  width-to-boundary  layer  thickness  ratio,  — 

8 

was  5.6  for  the  first  strip  heater  and  maintained  at  3  for  the  second  and  third 
heaters. 


The  preliminary  experiments  indicated  that  for  a  fixed  input  frequency  and 
power  level,  the  first  heater,  s  =  9,S  cm,  was  far  more  effective  in  forcing  the 
free  shear  layer  than  the  second  and  third  heaters,  set  at  s  =12  cm  and  s  =  15,5 
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cm,  respectively.  Also,  it  was  noticed  that  having  the  same  ratio  -7-;  the 
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second  heater  operated  more  effectively  than  the  third  heater.  Figure  59  reveals 
that  (T-S)  waves  initiated  by  the  first  and  second  heaters  benefitted  from  a  posi¬ 
tive  pressure  gradient  as  well  as  a  longer  period  of  amplification  than  that  of  the 
third  heater. 

In  the  final  version  of  the  experimental  model,  a  single  strip  heater  was  flush 
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mounted  at  s  =  9.2  with  -7-  =  5.3. 

s 

5.3.  External  Forcing  of  the  Non-Oscillating  Cavity  Flow 

The  receptivity  of  the  cavity  shear  layer  to  externally  imposed  disturbances 
was  studied.  Two  different  cavity  widths  were  examined,  with  corresponding 
b/0  =  66<  6/®lmin  aJfid  6/9=  77  <  6/9i„is„.  Throughout  this  section  these 
selected  widths  will  be  called  case  A  and  B  respectively.  In  both  cases  edge  velo¬ 
city  was  fixed  at  C/,  =  22cm  Abc  .  Cavity  flow  response  was  measured  in  terms  of 
the  maximum  shear  stress  that  was  developed  in  the  shear  layer  at  a  station 
.254  cm  (.1  inch)  upstream  of  the  downstream  corner. 

Spectral  analysis  of  the  response  velocity  fluctuations  indicated  a  one  to  one 
relation  between  forcing  and  the  response  frequency.  For  Case  A,  figure  60 
presents  the  response  of  the  cavity  sheair  layer  to  a  wide  band  of  the  forcing  fre¬ 
quencies  at  different  levels  of  forcing  power.  The  results  indicate  that  the  level 
of  sheeu*  layer  response  increases  as  the  forcing  power  increases.  For  forcing 
levels  above  9  watts,  resonance  peaks  appear  at  the  forcing  frequency  of  7.2  Hz. 


An  overall  phase  measurement  revealed  that  Case  A  satisfies  the  phase  cri¬ 
terion  -TT"'*'  2.  at  the  resonance  (F  =  7.2  Hz), thus  verifying  it  as  the 
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second  mode  of  oscillation.  The  absence  of  forcing  caused  the  resonance  mode 
to  disappear,  thus  revealing  that  the  flow  oscillations  are  not  self-sustaining. 
Figure  60  also  reveals  that  the  existence  of  resonance  depends  on  a  minimum 
amplitude  threshold. 

Figure  61  presents  the  response  of  the  cavity  shear  layer  for  Case  B  at  a  fixed 
forcing  level  of  9  watts.  This  case  shows  a  resonance  peak  at  the  forcing  fre¬ 
quency  of  6.2  Hz  with  2.  Appearance  of  resonance  peak  for  Case  B 
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at  a  lower  forcing  level  than  that  of  case  A  indicates  that  the  threshold  level 
decreases  as  b/do  increases.  In  Case  B,  the  effect  of  initial  lower  amplitude  of 
the  forcing  was  compensated  by  longer  cavity  length  in  such  a  way  that  reduced 
the  required  forcing  level  for  the  resonance  condition.  The  effect  of  higher  forc¬ 
ing  power  on  the  exponential  growth  of  the  perturbations  was  to  increase  the 
general  level  of  fluctuation  growth  (figure  62a).  Figure  62b  presents  the  growth 
behavior  of  the  different  frequencies  for  a  fixed  forcing  level  (9  watts)  for  Case 
B. 

It  can  be  concluded  that  external  frequencies  can  be  imposed  on  naturally 
non-oscillating  cavity  flows,  within  the  receptivity  range  of  the  shear  layer. 
These  frequencies  can  cause  the  cavity  flow  to  resonate  if  they  satisfy  the  follow¬ 
ing  conditions: 

1)  ^  for  the  given  width  (b):  and 

2)  A  >  Atjj  where  A  is  the  amplitude  of  external  frequency  at  the  doYm- 
stream  corner  and  Ath  is  the  threshold  amplitude.  It  is  logical  to  propose  that 
as  the  cavity  width  increases,  the  threshold  level  decreases  to  such  an  extent 


that  a  flow  backgroxmd  frequency,  which  satisfles  the  phase  criterion  and  has 
sufficient  amplitude,  will  initiate  the  self-sustained  oscillation. 


It  is  now  clear  that  due  to  the  absence  of  proper  frequencies  with  sufficient 
intensities  in  the  flow  background,  the  first  mode  did  not  appear  naturally  in 
this  study.  However,  by  selecting  a  proper  forcing  frequency  and  initial  ampli¬ 


tude,  the  first  mode  of  oscillations 


tt-®  1  simulated  for  Case  A. 

2ff 


This  was  accomplished  at  a  forcing  frequency  of  4  Hz.  a  power  level  of  35  watts. 
A  flow  visualization  of  this  case  is  shown  in  figure  63b,  In  comparison  to  figure 
63a.  which  is  the  same  case  but  unforced,  two  features  are  noteworthy.  The  first 
is  the  appearance  of  a  small  amplitude  wave  in  the  shear  layer  of  the  forced 
case.  The  second  is  the  interaction  of  this  wave  with  the  downstream  corner  and 
the  resultant  strong  recirculating  fiow  inside  the  cavity. 


It  is  interesting  to  note  that  the  response  of  the  cavity  flow  to  external  forc¬ 
ing  is  not  instantaneous.  Figures  64a  and  64b  present  simultaneous  traces  of 

the  velocity  fluctuations  taken  at  stations  corresponding  to  —  =  10  and  —  =  60. 

respectively.  The  cavity  was  in  the  non-oscillating  mode  with  6/9  =  66.  Figure 
62a  indicates  that  for  an  external  forcing  frequency  of  7  Hz.  simulating  Mode  11. 
it  takes  3  cycles  of  the  oscillation  before  the  peak  amplitude  establishes  itself. 
When  the  forcing  is  turned  off  it  takes  many  more  cycles  for  the  amplitude  of 
oscillations  to  diminish. 


5.3.1.  Effect  of  External  Forcing  on  Cavity  Drag  in  the  Nan-Oscillating 
Mode.  Case  B.  at  the  resonance,  had  a  drag  coefficient  of  0.0012.  This  drag 
coefficient  was  one  order  of  magnitude  higher  than  that  of  the  same  case  but  in 

the  absence  of  forcing.  Distributions  of  — g-  and  — ^  along  the  y  =  0  line  for 

c/,  t/, 

both  the  forced  and  unforced  cavity  are  given  in  figures  65a  and  65b.  These 
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figures  reveal  that  the  difference  in  the  drag  coefficient  of  forced  and  unforced 
cases  can  be  contributed  to  the  development  of  a  strong  shear  stress  level  in  the 
forced  case. 

5.4.  Effect  of  External  Fbrciiig  on  a  Naturally  Oscillating  ChvityFIow  ^b/9  >  b/&\^r, 
Once  the  cavity  flow  established  the  self-sustained  oscillations 
|6/i^  >  ft I  min]'  the  flow  resisted  any  external  disturbances.  Thus  in  order  to 

interact  with  the  flow  oscillations,  it  was  expected  that  the  amplitude  level  of 
external  forcing  must  be  comparable  to  that  of  the  self-sustained  fluctuations 
in  the  flow. 


A  non-dimensional  cavity  width  of  b/&o  =  95  was  selected  for  an  external 
forcing  experiment.  The  cavity  flow  was  oscillating  at  5.4  Hz  in  Mode  II.  An 
external  forcing  frequency  of  6  Hz  was  introduced  to  the  shear  layer.  Depicted 
in  figures  66a  through  66g  are  the  spectral  distributions  of  velocity  fluctuations 
for  several  levels  of  forcing  as  the  forcing  power  was  changed  from  2  to  27  watts. 

Spectra  reveal  that  as  the  amplitude  of  the  forcing  increases,  the  amplitude 
of  natural  oscillation  decreases  and  eventually  disappears.  At  the  forcing  power 
of  40  watts,  the  forced  frequency  becomes  the  dominant  frequency  of  the  oscil¬ 
lations.  Oscilloscope  traces  of  the  signal  for  the  case  when  both  natural  and 
forced  frequency  have  the  same  amplitude,  as  seen  in  figure  66d,  are  given  in 
figure  66h.  Strong  modulation  of  the  signal  indicates  that  both  frequencies  are 
present  simultaneously.  Note  that  in  the  co-existence  region,  figures  66b 
through  66f,  the  peak  amplitude  of  both  frequencies  is  lower  than  that  of  the 
single  frequency  regions  of  figures  66f  and  66g. 

Miksad  (1973)  in  his  study  of  dual  excitation  of  a  free  shear  layer  observed  a 
similar  interaction  of  two  frequencies.  He  reported  that  if  the  excitation 
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amplltudes  were  adjusted  so  that  both  disturbances  reached  finite  amplitude 
together,  the  final  amplitudes  of  both  components  tended  to  be  reduced,  but 
not  as  strongly  as  when  one  component  was  clearly  dominant.  Miksad  concluded 
that  the  distortion  of  the  mean  flow  by  the  stronger  frequency  made  velocity 
field  unsuitable  for  the  growth  of  the  weaker  frequency.  Similar  behavior  had 
also  been  observed  by  Sato  (1970)  in  S5niimetric  wakes. 

Important  to  our  study  is  the  fact  that  the  amplitude  of  the  forced  frequency 
in  figure  66g,  where  natural  frequency  Just  disappeeurs,  is  equal  to  that  of  the 
natural  frequency  in  the  absence  of  the  forced  frequency  (figure  66a).  This  fact 
establishes  the  amplitude  of  the  natural  frequency  as  the  level  which  should  be 
exceeded  by  the  amplitude  of  forced  frequency  to  change  the  frequency  of  the 
oscillations.  With  this  rule  it  is  possible  to  impose  forced  frequencies  on  the  cav¬ 
ity  sheaur  layer  within  the  receptivity  range  of  the  shear  layer.  These  forced 
oscillations  are  not,  in  general,  self-sustaining.  During  the  mode  switching 
region  (hysteresis  region),  self-sustaining  oscillation  in  a  desired  mode  could  be 
achieved  by  providing  a  proper  frequency  within  the  hysteresis  region  of  that 
mode.  Thus,  the  mode  switching  phenomena  could  be  eliminated  by  external 
forcing.  Transition  to  the  wake  mode  could  be  prevented  by  enhancing  the 
amplitude  of  natural  frequency  through  external  forcing.  It  is  possible  to 
extend  Mode  III,  last  mode  in  this  study, by  a  20  momentum  thickness  through 
external  forcing.  However,  once  the  cavity  entered  the  wake  mode,  it  was  not 
possible  to  reestablish  cavity  flow  oscillations  by  external  forcing. 

A  map  of  possible  interaction  with  the  cavity  flow  oscillations  observed  in  this 
experiment  is  given  in  figure  67. 
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5.5.  Cancdlatioii  Experiment 

Cavity  vridth  of  =  82,  with  frequency  of  6  Hz.  was  selected  for  the  cancel¬ 
lation  experiment.  The  signal  from  a  downstream  location  was  fed  to  a  function 
generator  (Tektronix  FG504)  which  had  a  built-in  phase  lock  loop.  It  was  possi¬ 
ble  to  lock  to  lower  harmonics  of  the  input  signal.  Phase  variation  of  the  output 
signal  with  respect  to  the  first  subharmonic  of  the  input  signal  was  possible 
through  an  independent  manual  adjustment  knob.  The  power  amplifier  was  then 
supplied  by  the  output  of  the  function  generator.  The  experimental  set  up  is 
sketched  in  figure  68. 

The  following  procedure  was  conducted  in  the  cancellation  experiment. 

1)  Output  amplitude  of  function  generator  was  set  at  0. 

2)  Phase  lock  loop  was  locked  to  first  subharmonic  of  the  input  signal . 

3)  Knowing  the  threshold  amplitude,  the  power  was  fixed  at  a  level 
which  was  enough  to  simulate  conditions  as  depicted  in  figure  66d. 

4)  Step  (3)  was  repeated  for  different  phase  angles  between  input  and 
output  of  the  function  generator  (phase  lock  loop). 

Due  to  the  small  rapid  phase  changes  of  the  cavity  flow  oscillations  it  was  not 
possible  to  have  a  phase  lock  situation  for  an  extended  period  of  time.  Occa¬ 
sionally  a  correct  phase  angle  could  be  found  where  the  forced  and  natural  fre¬ 
quencies  would  interact  in  such  a  way  that  wide  band  rms  fluctuations  were 


reduced.  A  sample  of  such  a  cancellation  is  shown  in  figure  69.  A  reduction  in 
rms  fiuctuation  was  achieved  by  a  factor  of  2  during  30  cycles  of  the  cavity 
shear  layer  oscillations,  liepmann  and  Nosenchuck  (19B2)  report  a  factor  of  2-3 
in  the  reduction  of  rms  fluctuations  of  the  natural  (T-S)  waves  in  the  fiat  plate 
boundary  layer. 


Chapter  6 


CONCLUSIONS 


6.1.  Summary  ct  Resulta 

(a)  Investigation  of  the  flow  fleld  in  the  oscillating  mode  reveals  that  the 
momentum  thickness  of  the  cavity  shear  layer  grows  linearly  with  x/Bo-  The 
linear  growth  of  the  momentum  thickness  may  indicate  the  existence  of  self¬ 
similarity  in  the  flow.  However,  the  self-similarity  argument  for  the  free  shear 
layer  requires  a  constant  maximum  shear  stress  along  the  layer.  In  contrast, 
cavity  shear  layers  show  an  exponentieJ  growth  of  the  maximum  shear  stress 
with  x/Bfj,  The  cavity  shear  layer  in  the  non-oscillating  mode  does  not  develop 
any  noticeable  shear  stress  along  its  span.  The  velocity  measurements  also  ver¬ 
ify  the  existence  of  a  recirculating  flow  inside  the  cavity.  Streamwise  and 
transverse  velocity  fluctuations  as  high  as  those  of  the  shear  layer  are  detected 
in  the  recirculating  region. 

(b)  Upstream  laminar  boundary  layer  has  a  shape  factor  value  of  2.5.  The 
effect  of  the  cavity  on  this  profile  is  found  to  decrease  the  shape  factor  to  a 
value  of  1.54,  which  is  close  to  that  of  the  turbulent  boundary  layer. 

(c)  The  drag  measurement  results  presented  in  figure  39  indicate  that  the 
drag  coefllcient  of  an  axisymmetric  cavity  with  laminar  upstream  boundary 
layer  remains  low  as  long  as  the  cavity  flow  stays  in  a  non-oscillating  or  self- 
sustained  oscillating  mode.  Results  also  reveal  that  the  drag  of  the  cavity  in  the 
non-oscillating  region  is  less  than  the  case  in  which  the  cavity  was  replaced  by 
a  solid  surface. 
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(d)  Flow  visuedization  study  (figure  18)  in  conjunction  with  the  drag  measure¬ 
ments  reveals  that  for  a  long  enough  cavity  width,  b/Ii>  1.2,  cavity  flow  exhibits 
a  wake-t3T>e  behavior. Also,  the  average  drag  coefficient  in  this  mode,  Cp  «  0.35, 
is  one  order  of  magnitude  higher  than  the  highest  drag  coefficient  in  the  oscil¬ 
lating  mode,  Cj}  =  .036. 

(e)  The  most  important  result  is  presented  in  figure  40.  The  drag  coefficient 
of  the  cavity  exhibits  an  exponential  dependence  with  6 /Ob  in  Mode  II  of  the 
oscillations  followed  by  a  decrease  in  the  same  mode  and  another  course  of 
exponential  growth  in  Mode  III.  An  independent  estimate  of  the  drag  coefficient 
through  measurement  of  the  shear  stress  and  mean  momentum  transfer  terms 
uV  along  the  y  =  0  line  confirmed  the  exponential  growth  of  the  cavity  drag. 
Furthermore,  these  measurements  revealed  that  the  uv  term  is  responsible  for 
30%  to  80%  of  the  net  dreig  coefficient  of  the  cavity. 

(f)  Study  of  the  cavity  pressvire  distribution  revealed  that  the  pressure  rise 
close  to  the  outer  region  of  the  downstream  comer  is  responsible  for  the  drag 
coefficient  of  the  cavity. 

(g)  In  the  oscillating  modes  of  the  cavity  flow,  the  maximum  pressure 

coefficient  (Cp„„)  always  occurs  at  the  upstream  corner  and  sta)^  so  until  the 
wake  mode  appears.  In  the  wake  mode  the  position  of  maximum  pressure 
coefficient  moves  toward  the  center  of  the  cavity  as  b/eo  increases. 

(h)  The  existence  of  the  feedback  loop  introduces  an  extra  stability  to  the 

cavity  flow  field  which  includes  a  fixation  feature  of  at  the  upstream 

comer  and  a  resultant  low  drag  for  the  cavity. 

(i)  The  distribution  of  the  shear  stress  along  a  straight  line  connecting  two 
comers  of  the  cavity  is  shoYm  to  have  direct  effect  on  the  positive  transverse 


velocity  (-H;)  distribution  on  the  same  line.  It  is  also  sho-wn  that  the  deflection 
of  the  shear  into  the  cavity  is  related  to  the  distribution  of  positive  transverse 
velocity  through  the  mass  balance  of  the  cavity. 

(j)  Stronger  deflection  of  the  shear  layer  into  the  cavity  causes  a  higher  ilv 
region  to  develop  close  to  the  dovmstream  corner  and  the  dreig  of  the  cavity  to 
increase.  Thus,  shear  stress  development  along  the  shear  layer  has  a  direct 
effect  on  the  drag  of  the  cavity  as  well  as  an  indirect  effect  through  its  role  on 
the  development  of  the  uv  term. 

The  type  of  shear  layer  interaction  with  the  dovnistream  comer  which 
depends  on  the  shear  layer  deflection  is  decisive  in  the  development  of  the  shear 
stress  throughout  the  cavity  shear  layer.  Based  on  the  findings  of  previous 
experiments,  it  is  proposed  that  inward  displacements  of  the  shear  layer  with 
respect  to  the  doYmstream  corner  will  result  in  an  overall  drop  of  the  shear 
stress  level.  As  a  result  of  this  overall  shear  stress  drop,  a  less  positive  region  of 
transverse  velocity  develops.  Consequently,  the  mass  balance  of  the  cavity 
requires  less  deflection  of  the  shear  into  the  cavity.  It  also  requires  that  the 
shear  layer  deflect  back  and  restore  the  position  of  the  shear  layer.  This 
mechanism  keeps  the  position  of  fixed  at  the  upstream  corner  and  it 

might  also  introduce  some  low  frequency  components  to  the  oscillations. 

(k)  The  control  of  the  cavity  flow  oscillation  is  facilitated  by  external  forcing 
of  the  cavity  shear  layer.  It  is  shown  that  the  excited  T-S  waves  can  be  used  to 
force  the  shear  layer  generated  from  separation  of  respective  boundary  layers. 
Thus,  a  transformation  of  the  T-S  waves  of  the  boundary  layer  to  the  Kelvin- 
Helmholtz  waves  of  the  shear  layer  is  possible. 

(l)  For  short  cavities,  i.e.,  non-oscillating  mode  6 /9o  <  82,  periodic  flow  oscil¬ 
lation  at  a  wide  range  of  frequencies  was  obtained.  A  resonance  condition  was 


observed  for  a  selected  forcing  frequency  if  it 

1)  satified  the  phase  criterion  N 

ecid 

2)  its  amplitude  exceeded  a  certain  threshold  amplitude. 

For  a  fixed  amplitude  the  threshold  level  increases  as  b  decreases.  The  forced 
resonance  cases  are  not  self-sustaining. 

(m)  Previous  conclusions  indicate  that,  for  a  given  6/(9o.  the  first  frequency 
in  the  fiow  background  that  satisfies  the  phase  criterion  and  has  sufficiently 
high  amplitude  will  start  the  self -sustained  oscillations.  Thus  the  determination 
of  the  h/dolmin  strongly  depends  on  the  turbulence  content  of  the  flow  back¬ 
ground. 

(n)  Cavity  drag  in  the  non-oscillating  mode  can  be  increased  by  one  order  of 
magnitude  through  external  forcing.  Non-oscillating  cavity  flow  at  6/9o  =  77 
produces  a  drag  coefficient  value  of  1.3  x  10~^.  External  forcing  at  6  Hz  and 
power  level  of  9  watts  increases  the  drag  coefficient  to  a  value  of  1.2  x  10“®.  It  is 
also  possible  to  generate  shear  stress  levels  several  orders  of  magnitude  higher 
than  that  of  the  laminar  shear  layer  by  external  forcing. 

(o)  For  naturally  occurring  oscillations  it  is  possible  for  two  waves  to  co-exist 
in  the  shear  layer  (natural  and  forced).  Spectral  analysis,  however,  shows  that 
the  amplitude  of  the  forcing  frequency  has  a  direct  effect  on  the  growth  of  the 
natural  frequency.  External  forcing  made  possible  the  ability  to  completely  sub¬ 
sidize  the  natural  self-sustained  frequency  by  a  non-self-sustained  frequency.  It 
is  shown  that  the  amplitude  of  the  established  self-sustained  oscillation  is  the 
amplitude  levc’  that  should  be  exceeded  by  the  amplitude  of  the  forced  fre¬ 
quency  for  complete  switching  of  self-sustained  to  forced  frequency. 


(p)  Mode  switching  or  transition  to  the  wadce  mode  can  be  delayed  by  enhanc¬ 
ing  the  self-sustained  oscillations  through  external  forcing.  Once  the  weike 
mode  appears,  the  self -sustained  oscillations  cannot  be  reestablished.  Different 
mode  extension  and  simulation  is  possible  by  providing  the  proper  forcing  fre¬ 
quency  to  the  cavity  shear  layer. 

(q)  The  most  important  result  of  external  forcing  is  the  possibility  of  cancel¬ 
ling  or  damping  the  established  self-sustained  oscillations.  An  external  forcing 
frequency  equivalent  to  that  of  the  cavity  with  an  amplitude  less  than  threshold 
level  which  has  a  proper  phase  angle  respective  to  that  of  the  self-sustained  one 
was  able  to  reduce  rms  fluctuations  by  a  factor  of  2. 

6.2.  Concluding  Remarks 

This  study  has  disclosed  that  the  presence  of  self-sustained  oscillations  helps 
in  keeping  a  low  level  of  drag  coefficient  for  the  cavity.  A  more  detailed  study  of 
the  role  of  the  downstream  corner  on  the  overall  development  of  the  shear 
stress  along  the  shear  layer  and  on  the  drag  of  the  cavity  is  desirable.  It  is 
important  to  clarify  the  effect  of  the  cavity  presence  on  the  net  drag  coefficient 
of  the  system. 

With  respect  to  the  active  control  of  the  cavity  flow  oscillations,  if  a  more  ela¬ 
borate  feed-back  system  is  used,  a  high  percentage  of  fluctuation  cancellation 
might  be  attainable. 


Appendix  A 


Determination  of  Overall  Riase  Difference  for  the  Cavi^ 

For  the  purpose  of  obtaining  phase  information  between  two  time  dependent 
signeds,  the  cross-correlation  function  between  these  two  signals  should  be  car¬ 
ried  out.  By  cross-correlating  two  signals,  the  relative  non-dimensional  phase 
change  between  the  location  of  the  reference  probe  and  the  second  probe  can  be 
determined.  This  is  done,  as  shown  in  sketch  (A-1),  by  measuring  the  time 
shift,  along  the  positive  r-eucis,  between  two  corresponding  peaks  in  the  cross- 
correlation  function  of  two  signals  and  dividing  this  time  by  the  average  period 
between  peaks. 


SKETCH  A-1 


Then  the  non-dimensional  phase  difference  between  two  locations  separated  by 
Ax  will  be 


StT  TJ^ygJ{J^Qg  PERIOD 


The  wavelength  of  the  oscillation  can  be  determined  as 


X  =  Ax 


grr 

A^ 


(2 


Knowing  the  fundamental  frequency  of  oscillation,  the  phase  velocity  can  be  cal 
culated  from 


£/e  =/X 


(3 


Overall  phase  difference  between  the  two  corners  is  given  by 


-66- 


Appendix  B 

Laser  Dof^ler  YdocLzneter  Optical  Arraneement 

A  two  watt  argon-ion  laser  Is  used  as  the  light  source.  Past  a  collimating  lens, 
a  color  separator  provides  two  lines  of  the  laser,  the  514  nm  green  line  and  488 
nm  blue  line.  Each  beam  is  split  into  two  pedrs  of  equal  intensity  beams  by  a 
polarization  rotator  and  a  beam  splitter.  Each  pair  passes  through  its  desig¬ 
nated  frequency  shifter.  Non-desired  diffracted  beams  are  blocked  by  two  slit 
mechanisms.  The  measuring  volume  is  formed  by  focussing  the  four  parallel 
beams  by  a  single  598  mm  focal  length  lens.  The  collimator  is  used  to  minimize 
the  lateral  extent  of  the  focal  volume. 

The  scattered  light  from  the  particles  is  collected  by  a  600  mm  focal  length 
lens.  A  dichroic  mirror  is  used  for  color  separation,  with  a  focussing  lens  used 
to  focus  the  scattered  light  of  each  of  the  two  colors  into  a  pin-hole  (350  jj,  m) 
and  photomultiplier  tube  assembly.  A  first-surface  flat  mirror  is  used  to  rotate 
the  beam  90“,  A  second  mirror  is  used  to  deflect  the  scattered  light  into  the  col¬ 
lecting  optics. 

The  component  of  the  velocity  vector  measured  by  the  LDV  and  the  Doppler 
constant,  the  conversion  factor  from  frequency  to  velocity,  is  determined  by  the 
transmitting  optics.  The  green  beams' pair  enters  the  test  section  in  a  plane 
forming  an  angle  of  0®  with  the  downstream  direction,  while  the  blue  beauns' 
plane  forms  an  angle  of  90®.  Therefore,  streamwise  velocity  component  is  meas¬ 
ured  by  the  fringe  system  formed  by  the  green  beams  while  velocity  components 
perpendicular  to  the  mean  free  stream  are  measured  by  the  fringe  system 
formed  by  the  blue  beams.  Figure  B-1  presents  the  schematic  of  focal  volume 
formation.  The  Doppler  constant  which  is  related  to  the  fringe  spacing  within 


the  focal  volume  is  given  by  the  equation 


K  = 


X 


2  sin 


2 


where  X  is  the  light  wavelength  and  a  the  angle  at  which  the  beams  cross  in  the 
focal  volume.  The  angle  a  =  4.72“  was  used  which  gives  Doppler  constants  for 
the  green  and  blue  components  kq  =  6.24r  {m^ec)  (Si Hz  )~^  and 
Kb  =  5.92  (m^ec)  {liHz)~^  respectively. 


Appendix  C 


Velocity  Bias  and  Bragg  Cell 

Fringes  of  the  measuring  volume  are  formed  by  the  interfering  effect  of  the 
two  beams  crossing.  These  fringes  are  stationary  and  do  not  produce  Doppler 
frequency  for  a  stationary  particle  in  the  measuring  volume. 

Shifting  the  frequency  of  one  beam  can  produce  a  system  of  moving  fringes. 
Then,  a  stationary  particle  in  the  measuring  volume  will  give  an  output  fre¬ 
quency  from  the  photodetector.  If  the  velocity  increases  in  the  direction  of  the 
fringe  movement,  the  frequency  will  decrease.  If  it  is  opposite  to  the  fringe 
movement,  the  frequency  will  increase.  Therefore  knowing  the  direction  of  the 
fringe  system  movement  will  discern  the  particle  direction. 

In  this  experiment,  two  acousto-optic  cells  (Bragg  cell)  were  used  to  fre¬ 
quency  shift  one  beam  of  each  color  by  a  fixed  40  MHz.  Each  cell  is  installed  in  a 
cylindrical  housing  that  aligns  with  the  optical  axis.  A  compensating  wedge  is 
used  to  correct  for  shifted  beam  deflection,  and  a  tilt  adjustment  allows  shifting 
light  intensity  to  the  deflected  beam,  as  much  as  80%.  The  cell  is  derived  by  a  40 
MHz.  2  watt  signal  provided  by  a  power  amplifier.  A  frequency  of  40  MHz  is  often 
too  large.  A  downmix  circuit  permits  decrease  of  the  effective  frequency  shift  of 
the  photodetector  output  to  values  from  2kHz  to  10  MHz  in  a  2,  5,  10.  sequence. 
Breigg  shifts  of  nominally  30  and  50  kHz  were  used.  In  this  study  the  green  (U 
component)  channel  was  shifted  against  the  mean  flow,  while  the  blue  channel 
(V  component)  was  shifted  in  the  positive  y  direction  (figure  B-1). 
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Appendix  D 


IJ}V  Signal  Processiiig 

The  burst  processing  encompasses  the  required  electronics  for  the  actual 
measurement  of  the  particle  velocity.  A  particle  burst  is  considered  valid  if  it 
satisfies  the  following  conditions; 

1)  Its  amplitude  crosses  zero  and  exceeds  the  threshold  level  (about  50  mV) 

2)  For  zero  crossing  (with  h5^teresis),  the  signal  should  be  more  than  10  mV. 


The  pulse  that  occurs  on  a  zero  crossing  after  the  threshold  has  been  exceeded 
is  the  signal  of  the  interest.  At  this  point  the  leading  edge  of  a  time  envelope  is 
marked  (by  the  processor)  and  gates  a  high  speed  clock  (250  MHz)  into  a 
counter.  The  length  of  the  envelope  is  determined  by  the  total  number  of  cycles 
exceeding  the  threshold.  The  first  cycle  that  exceeds  the  hysteresis  level  but 
not  the  threshold  level  determines  the  trailing  edge  of  the  envelope.  The  trailing 
edge  stops  the  counter  and  the  clock.  The  input  frequency  is  then: 


INPUT  FREQUENCY  (a) 


ENVELOPE  CYCLE  NUMBERS 
TOTAL  ENVELOPE  TIME 


The  resolution  of  time  measurement  is  1  nsec.  In  this  experiment  only  bursts 
with  more  than  16  cycles  were  used  for  processing.  A  schematic  of  valid  burst 
processing  is  given  in  figure  D-L  The  output  of  the  burst  processor  is  in  digital 
form.  It  feeds  into  the  computer  interface  that  incorporates  the  coincidence 
test  between  the  two  processors.  The  coincidence  test  consists  of  the  initiation 
of  a  time  window  after  completion  of  processing  by  one  of  the  processors.  If 


data  from  the  other  processor  are  not  received  before  the  end  of  the  window,  the 
data  are  discharged;  otherwise  coincidence  is  attained  and  the  data  from  both 
processors  are  transferred  to  the  computer. 


A  minicomputer  was  used  to  acquire  and  store  the  data.  A  computer  program 
was  developed  for  data  acquisition  using  the  LDV  S5rstem.  This  computer  pro¬ 
gram  accomplishes  two  tasks.  It  generates  a  header  which  contains  information 
relevant  to  the  conditions  of  the  experiment  and  the  LDV  setup.  It  also  acquires 
a  buffer  of  16,000  words  of  data.  This  buffer  size  corresponds  to  a  total  of  3200 
velocity  measurements  in  two-component  measurements,  since  two  words  are 
generated  for  each  component  and  one  word  was  used  for  the  time  interval 
information  between  data  points. 


Appendix  E 


UJVData  Processing 

The  following  expressions  are  used  to  calculate  the  mean,  second  moment 
and  cross-correlation  between  two  velocity  components. 
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where  At*  is  the  time  between  data  points  and  cj  is  the  frequency.  The  sub¬ 
scripts  G  and  B  refer  to  the  frequencies  measured  by  the  green  and  blue 
beams,  respectively.  The  equations  relating  the  mean  and  the  correlations  of  the 
various  velocity  components  are 

U  =  |c5c  -  /’5c]  Kg,  V  =^g  -  /’5jj|  Kg  (E-6) 

=  [z?c  -  FSc  +  /’52c]  (*^-7) 

^  =  \l?g  -  2ag  FSg  +  ^51]  K^g 

wS  =  ^Sc®S  ~  OoFSg  —  OgFSc  +  /’5c/’5j5]  Kg  X  Kg  (E-9) 


where  FSq  and  FSg  are  frequency  shifts  of  green  and  blue  channels, 
respectively,  and  k  is  the  frequency  to  velocity  conversion  factor. 

Rms  fluctuations  and  shear  stress  values  can  be  obtained  using  the  following 
relationships: 


CT  = 
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Fig.  6  Relative  location  of  the  LDV  with  respect  to  the  model  and  characteristics  of  the 
measuring  volume  (Table  2). 
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stream  velocity. 


(a)  Mode  II,  F  =  6  Hz,  vertical  scale  .02  volt/div 


(b)  Transition  F  =  4.6  Hz  Mode  II 

F  =  6.4  Hz  Mode  III,  Vertical  scale  .02  volt/div 


(c)  Mode  111  F  =  6  Hz  vertical  scale  .04  volt/div 


Oscilloscope  traces  of  cavity  shear  layer  oscillations 
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Fig.  19  Non-dimensional  mean  streamwise  velocity  profiles 
F  =  6  Hz,  Mode  III. 
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)  where  u/U  -  0.5  for  cavity  widths  q  b/G[.  =  85  and 


Fig.  30  Growth  of  the  maximum  shear  stress  in  the  cavity  shear  layer 
with  downstream  distance; 
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nd  downstream  face,  respect ively. 


Fig.  48  Continued 
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Fig.  52  Distribution  of  the 
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Fig.  55  Distribution  of  shear  stress  along  y  *  0  line  in  Mode  III 
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Fig.  66  Velocity  fluctuation  spectra,  presenting  interaction  of  the 
natural  and  forced  frequencies  in  the  cavity  shear  layer  for 
different  power  levels  (P,),  (h)  scope  trace  of  Case  (d) . 
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Fig.  67  Possible  interactions  with  the  cavity  flow  through 
external  forcing 

1)  Mapping  of  Mode  II  into  Regions  1  and  2. 

2)  Mapping  of  Mode  III  into  Region  3. 

3)  Extension  of  Modes  II  and  III. 

4)  Delay  of  the  wake  Mode. 

5)  Simulation  of  Mode  I. 
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Fig.  68  Schematic  of  the  cancellation  set-up. 


Fig.  69  Damping  of  natural  oscillations  through 
external  forcing 

(a)  natural  oscillation 

(b)  strip-heater  on. 


